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1. Introduction 

Polyelectrolytes are polymers bearing ionizable groups, which, in polar solvents, can 
dissociate into charged polymer chains (macro-ions) and small " counter-ions" . Well known 
examples of such systems are proteins, nucleic acids or synthetic systems such as sulfonated 
polystyrene or polyacrylicacid. Very often, the polar solvent is water. Because of their 
fundamental importance in biology and biochemistry, and because of their hydrosolubility, 
ionizable polymers have been the object of a continued interest since the early days of 
polymer science (see e.g. [1,2] and references therein). It is however, a widely acknowl- 
edged fact [3] that they remain among the least understood systems in macromolecular 
science. This situation is in sharp contrast with the case of neutral polymers. In fact, 
even the properties that are now routinely used to characterize neutral polymers (e.g. 
light scattering, osmometry, viscosity..) are still very poorly understood in polyelectrolyte 
solutions. 

Briefly speaking, the origin of this difference with neutral polymers can be traced back 
to the difficulty to apply renormalization group theories and scaling ideas to systems in 
which long range (coulomb) forces are present. The success of the modern theories (and, 
although this was not always realized as they were developped, of many older approaches) 
of neutral polymer solutions is based on the fact that the range of the interactions be- 
tween molecules is much smaller than the scale determining the physical properties of the 
solution, which is the size of the polymer chain or the correlation length. If, as it is in 
general the case in polymer physics, the main issue is the variation of solution properties 
with molecular weight, the interactions only affect prefactors that can be adjusted to ex- 
perimental results. The theories [4,5,6] developped in this spirit have proven extremely 
powerful for the interpretation of experimental data on neutral polymer solutions. Poly- 
electrolyte solutions are more complex, both short range (excluded volume) and long range 
(coulombic) interactions are simultaneously present. The screening of the electrostatic in- 

3 



teractions introduces an intermediate length scale in the problem, that can be comparable 
to the chain size or to the correlation length. Moreover, the details of the local chain 
structure are important and control the phenomenon of counterion condensation [7] (see 
section 3). This condensation in turn modifies the long range part of the interaction. The 
long range interactions actually also has have a nontrivial influence on the local structure 
[8,9,10] (stiffness) of the polymer chains (section 4). The implication of this complicated 
coupling between small and large length scales is that the theoretical results for polyelec- 
trolytes can be expected to be more "model-dependent" than for neutral polymers. In 
particular, the functional dependence of several physical properties on molecular weight or 
concentration depends in some cases on the local description that is used for the polymer 
chain. Such situations never occur in neutral polymers. 

Intrinsically related to this theoretical problem, is the more practical difficulty in com- 
paring experiments and theory for polyelectrolyte solutions. In neutral polymers, the 
comparison is made possible by adjusting interaction dependent prefactors and fitting ex- 
perimental results to the theoretically predicted dependence on, e.g. molecular weight. In 
polyelectrolytes, "microscopic" parameters can have bigger influence, and a much more 
precise modelling might be necessary in order to interpret experimental results. This of- 
ten involves the introduction of more adjustable parameters in the interpretation, which 
sometimes becomes inconclusive. 

In spite of these considerations, the modern approach to neutral polymers physics has 
inspired a number of contributions to the polyelectrolyte problem in the last two decades. 
By "modern", we mean here those theoretical approaches that attempt to use as much 
as possible "coarse grained" models, with a minimal number of microscopic parameters. 
The goal of this review is to summarize these approaches, with the hope that they could 
provide a useful conceptual framework for further developpment. The review is organized 
as follows. Sections 2 to 4 are concerned with single chain properties. Section 2 deals with 
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the academic problem of a charged chain in the absence of counterions. The phenomenon 
of screening, and the influence of the small ions on the interaction between monomers, 
is discussed in section 3. The influence of long range interaction on the local stiffness 
(electrostatic rigidity) is described in section 4. In sections 5 and 6, the properties of 
interacting polyelectrolyte chains are considered. The simpler case of gels and brushes 
-in which the structure is controlled by the preparation of the system- is described first. 
Semidilute and concentrated solutions are the subject of section 6. Finally, some results 
on dynamical properties (diffusion, electrophoretic mobility and viscosity) are presented 
in section 7. 

Some important subjects have been left out of the scope of this review, either be- 
cause they have been reviewed recently elsewhere or because our understanding is too 
fragmentary. Among those are the dielectric properties of polyelectrolyte solutions [11], 
their rheological properties, computer simulation results [12], and all the problems related 
to polyelectrolyte adsorption [13]. Finally, we systematically refer the reader to the origi- 
nal papers for detailed comparison with the experimental results. The reason, mentioned 
above, is that such comparisons are, in general, not simple, and involve the determination 
of several parameters. 

2. Charged chains at infinite dilution - asymptotic properties 

2.1 Definition of the Model and Flory-like calculation 

In the theoretical limit of infinite dilution and zero salt concentration, a polyelectrolyte 
chain can be modelled as a connected sequence of charged and uncharged monomers in a 
"dielectric vacuum" that represents the solvent. This model is not very realistic to describe 
actual polyelectrolyte solutions, as will become clear in our discussion of ion condensation 
and screening but it serves here to introduce some important concepts relative to charged 
polymer chains. 
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The energy for such an isolated polymer chain, made up of N monomers, is written as: 

where Hq is the energy of the neutral polymer chain that would be obtained by " switching 
off' the electrostatic interaction, and the double summation on the left hand side includes 
all monomer pairs. Here Z{ is the charge (in units of the electronic charge e) of monomer 
i, T is the absolute temperature and i B = e 2 /(47reoe r /c5T) is the Bjerrum length, that 
characterizes the strength of electrostatic interactions in the solvent. For water at T = 
300X, l B ~ 0.7nm. 

The "neutral polymer" energy Hq contains terms describing the chemical bonding of 
the monomers as well as short range excluded volume terms. In this section, we focus on 
the very simple case where Hq is the free energy of a gaussian, or bead-spring, chain, i.e. 

= ^2~ ^ ( r *- r *+l) • (2.2) 

i=l,JV-l 

The "monomer size" b is such that the mean squared averaged end-to-end distance for 
the neutral chain is Rq = Nb 2 . The physical situation described by this model is that of 
a very weakly charged, flexible chain. In that case, each "spring" actually represents a 
flexible subchain of several neutral monomers. Moreover, the absence of excluded volume 
interactions implies theta conditions for the neutral backbone. Although these conditions 
are rather unlikely to be achieved in experiments, the model contains the two key in- 
gredients that monitor the asymptotic behavior of charged polymers in the absence of 
screening, namely long range interactions and monomer connectivity. These are the only 
relevant contributions as can be seen from the form of the Flory free energy for a chain 
of N monomers, carrying a total charge fNe (i.e. / is the fraction of charged monomers, 
(1 — f)N monomers are neutral). The Flory energy [14] for a chain 1 of size R is the sum 



1 The Flory free energy ignores numerical prefactors of order unity; the prefactor of the 
electrostatic energy depends on the actual distribution of the charges in the sphere of 
radius R, it is equal to 3/5 if the charges are uniformly distributed 
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of an elastic energy kgTR 2 /Nb 2 and of a Coulombic energy kgT(N f) 2 £ B / R, 

E Flory ~ k B T + ^ J , (2-3) 

which upon minimization with respect to R gives an equilibrium size 

R~Nf 2 l\l B tf)W . (2.4) 

This establishes the well known result, that the chain size is proportional to the number 
of monomers and to a 1/3 power of the charge fraction /. This result was first obtained 
by Katchalsky and coworkers [15] using a distribution function approach, shortly before 
Flory published an equivalent calculation for chains with short range excluded volume. 

The electrostatic interactions tend to swell the chain, therefore the equilibrium radius 
of a weakly charged polyelectrolyte chain must be larger than the gaussian radius Rq. 
This condition defines a minimum charge fraction Jq ~ iV~ 3 / 4 (&//&) 1//2 above which the 
electrostatic interactions are relevant and the radius is given by (2.4). If the fraction of 
charges / is smaller than /q the polyelectrolyte chain essentially behaves as a neutral chain 
and has a gaussian radius Rq. In the following we always implicitly assume that f > fc- 

In contrast to Flory's calculation, whose success is known [5] to arise from an uncon- 
trolled cancellation of errors, Katchalsky's approach is now believed to yield an essentially 
exact scaling result, as will become clear in the next sections. A more refined calculation 
of the electrostatic energy, accounting for the fact that the shape of the chain is rodlike 
rather than spherical, only yields a logarithmic correction to (2.4) 

i? ~ iV/ 2 / 3 (^ B 6 2 ln(JV)) 1 / 3 . (2.5) 

Because of this rodlike character, the polyelectrolyte chain is often described as being a 
"fully extended" object. This statement, however, somewhat overinterprets the actual 
meaning of relation (2.4). Short range molecular flexibility does not have to be frozen in 
a locally stretched configuration to obtain such relations. 
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To understand this point, and to get a deeper insight into the meaning of Katchalsky's 
calculation, it is useful to introduce the concept of electrostatic blob [14]. An electrostatic 
blob is a chain subunit within which the electrostatic interactions can be considered as a 
weak perturbation. The size £ e and the number of monomers g in such a subunit are thus 
related by 

Uaf^B _ 1 (2.6) 



He 

Another relation between £ e and g is obtained from the gaussian statistics of the chains 

on the scale £ e , £g ~ gb 2 . This gives 

£e ~ b(f X/b)- 1 / 3 

(2.7) 

The size of the chain (equation (2.4)) is therefore (N/g)^ e , so that the chain can be thought 
of as a linear (in the sense that its lateral dimensions are vanishingly small compared to 
its length) string of "electrostatic blobs" (figure 1). 

The blob picture is particularly useful to discuss the structure of a charged chain in 
a good or poor solvent. In a good solvent, the only modification compared to the theta 
solvent case is the good solvent statistics characterized by a swelling exponent v ^ 1/2 
must be used in the relation that relates the blob size to the number of monomers in a 
blob, £ e ~ g h ' . The result is still a linear string of blobs, R ~ N f y , but the exponent y 
that describes the variation with the charge fraction is now (2 — 2v)/{2 — v). In a poor 
solvent, the problem is more subtle, and has been analysed in detail by Khokhlov [16]. In 
that case, the chain would be completely collapsed into a globular state in the absence 
of electrostatic interactions. The effect of these interactions is to break up the globular 
structure into a string of " globular blobs" . The criterion that determines the blob size is 
now that the electrostatic energy within the blob is sufficient to overcome the interfacial 
energy that must be paid when monomers are exposed to the solvent, i.e. 

(fg)X 



He 
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TZi/(k B T) (2.8) 



where F is the inter facial tension between the collapsed phase and the pure solvent. A 
second relation is obtained by expressing that the blob has the same structure as a small 
section of collapsed polymer, 

9 ~ <l . (2.9) 

Here n is the density of the collapsed chain. In the simple case where the second and 
third virial coefficients between monomers are, respectively, 6 3 (T — 9)/T = — 6 3 r and 6 6 , 
n ~ t&~ 3 and F = kj^TT^b~^. Equations (2.7) and (2.4) are replaced by 

~ b(f\/b)- 1/3 

g-rifX/b)- 1 (2.10) 
R-iN/g^e-Nr^iflJbfl* . 

2.2 Variational approaches 

The Flory-type calculation discussed above is a typical mean-field approach. Similar 
approaches which have been successfully applied to the charged bead spring model of 
section 2.1. include the "chain under tension" model [14,17] and the "gaussian variational 
method" [5,17,18,19]. Both methods are mean-field approaches, formulated in a variational 
fashion. The chain is approximated by a simpler system described by a trial Hamiltonian 
Ht, and the "best" trial Hamiltonian is obtained by minimizing the variational free energy 

F var =< H > t -TS t (2.11) 

where < H >t is the average value of the energy for the true system, taken with a statistical 
weight exp(— Ht/hgT), and St is the entropy of the trial system. 

In the 'chain under tension' model, the trial Hamiltonian is that of a noninteracting 
gaussian chain subject to an external tension F, 

Ht (F) = -^g- " r *) 2 + F -( r iV " ri) • (2.12) 

i=l 
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This model is simple enough to allow explicit calculations of many properties of the chain. 
In particular, it can be shown that the minimization of the corresponding variational free 
energy yields a force 



The results for the chain size and chain structure factor were shown to agree quantitatively 
with Monte-Carlo calculations [17]. This model can also be used to give a quantitative 
meaning to the "electrostatic blob size" introduced above. This is done by noting that 
in the chain under tension, the mean squared distance between two monomers separated 
by n monomers along the chain is i?(n) 2 = n& 2 + n 2 (6 2 F/3/c^T) 2 . The blob size can be 
obtained by identifying the second term with ((n/,g)£ e ) 2 , £ e = in agremeent with 

(2.7). The blob size is therefore directly related to the prefactor of n 2 in the expression of 



In the gaussian variational approach, the trial Hamiltonian is a general quadratic form 
of all monomers coordinates. This method, which is tractable for cyclic chains only, was 
first introduced by des Cloizeaux [5] as an alternative to Flory's theory for chains with 
short range excluded volume. In that case, however, the method gives a swelling exponent 
v = 2/3 in three dimensions, and is therefore inferior to Flory's approach. For long range 
interactions. It was shown [18] that the variational gaussian approach predicts the same 
asymptotic behaviour as Flory's theory and as the "chain under tension" method. As 
for the chain under tension model, the quantitative agreement with numerical simulation 
results is good [17]. 

2.3 Renormalisation group calculations 

The coincidence of the predictions from the Flory approach and the variational gaussian 




(2.13) 



and a chain size 




(2.14) 




10 



method suggests that the predicted asymptotic behaviour is indeed exact, except may be 
for logarithmic corrections. Exact renormalisation group calculations were performed for 
space dimensions d > 4 [20,21], showing that for 4 < d < 6 the swelling exponent of 
the chain is v = 2/(d — 2) (as does the variational calculation of reference [18]). No such 
calculation is available for d = 3, but the result for d = 4 suggests that the maximum value 
v = 1 is obtained for d < 4, while gaussian behaviour is recovered for d > 6. Moreover, 
a "naive" real space renormalisation argument [4] also predicts v = 1. Note that a naive 
extension of the Flory theory to d dimensions predicts an upper critical dimension where 
the chains are gaussian equal to 6 but a lower critical dimension where v = 1 equal to 3 
in disagreement with the renormalisation group results. 

The success of variational and mean-field methods for chains with long range inter- 
actions is not unexpected. As pointed out in [18], these methods usually fail because a 
poor approximation of the monomer-monomer correlation function is used in the calcula- 
tion of the potential energy. A gaussian correlation function will be maximum when two 
monomers are in contact, whereas the true correlation function vanishes in that case. For 
long range forces, most of the potential energy comes from configurations in which the 
monomers are far apart, so that the gaussian approximation becomes much better. 

2-4 Screening of electrostatic interactions 

When a finite concentration of salt is present in the solution - which is always the case 
in experiments, at least because of water dissociation-, the solution becomes a conducting, 
rather than dielectric, medium. The fundamental implication is that the Coulomb inter- 
action between charged monomers is screened by the salt solution, i.e. the electrostatic 
potential created by a monomer, or a group of monomers, falls off exponentially rather 
than algebraically with distance. Thus distant parts of the chain do not interact, and the 
chain can be expected to behave, in the asymptotic limit N — > 00, as a random walk with 
short-ranged repulsive interactions. 
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It is a general result from the theory of charged solutions that charge fluctuations 
become uncorrelated (or, equivalently, that an external charge distribution is screened) 
over a typical distance kJ^, where k s is related to the thermodynamic properties of the 
solution [22] . The screened interaction between monomers, mediated by the salt solution, 
can formally be obtained by integrating out the degrees of freedom associated to the 
salt ions in the partition function as shown in Appendix A. The resulting interaction is in 
general a complex function of the monomer coordinates, involving many body interactions. 
In the case of dilute salt solutions and sufficiently weak perturbations, however, the linear 
response (or Debye-Hiickel) theory can be used, and the salt can be treated as an ideal 
gas. The screening length k^ 1 is then given by 

k 2 = 4n£ B I, (2.15) 

where I is the ionic strength of the solution defined as I = Xlionic species ^i c i^ being 
the valence and q the concentration of species i. 

The effective interaction can in this limit be written as a sum of pairwise additive 
interactions between the monomers. The effective pair potential is given by the Debye- 
Hiickel formula [23] 

v DH (r) = k B T^exp(- K r) . (2.16) 
r 

This short range interaction potential is the starting point of many theoretical studies 
of polyelectrolyte solutions. The definition of an effective pair potential requires however 
two conditions. The first condition, that the ionic solution is dilute, is usually satisfied 
in polyelectrolyte solutions. The second condition of a weak external perturbation intro- 
duced by the polymer on the small ion solution must be considered more carefully, and is 
addressed in the next section. 

Although the potential (2.16) is short-ranged, its range k^ 1 can be much larger than 
the monomer size, which is the typical interaction range in neutral polymers. Depending 
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on the ionic strength, k can vary typically from less than lnm to more than lOOnm. 
At distances smaller than the screened potential (2.16) is close to the unscreened 

Coulomb potential discussed in section (2.1). The behaviour and properties of polyelec- 
trolyte chains can therefore be markedly different from those of neutral chains over rather 
large length scales. In fact, section 4 will show that in some cases the effect of electrostatic 
interactions extends over length scales much larger than k~ 1 . 

Finally, it should be noted that the concept of effective interaction between monomers, 
obtained by integrating out the degrees of freedom of the small ions, is useful only for 
the calculation of static properties. Dynamical properties such as diffusion constants or 
viscosities do not have to be identical to those of an hypothetical polymer solution in which 
the monomer-monomer interaction would be given by (2.16). The ionic degrees of freedom 
can play a nontrivial role in the dynamics, as they do in suspensions of charged colloids 
[24]. 

2.5 Annealed and quenched polyelectrolytes 

An important degree of freedom of polyelectrolyte chains is the distribution of the 
charges along the chemical sequence. This distribution can be either quenched or annealed 
[25]. When a weakly charged polyelectrolyte is obtained by copolymerisation of a neutral 
and of a charged monomer, the number of charges on each chain and the positions of 
the charges along the chain are fixed, the distribution is quenched. If the charges are 
sufficiently regularly distributed along the chain and do not have a tendency to form long 
blocks, the precise distribution of the charges does not seem to be a relevant variable for 
the chain statistics. At the level of scaling laws that we use, one can expect that only 
the numerical prefactors are different when polyelectrolytes with a random or a periodic 
distribution of charges are compared. 

Polyacids or polybases are polymers where the monomers can dissociate depending on 
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the pH of the solution and acquire a charge. The dissociation of an H + ion from a polyacid 
gives rise to the apparition of a COO group and thus of a negative charge. This is an 
annealed process, the total number of charges on a given chain is not fixed but the chemical 
potential of the H + ions and thus the chemical potential of the charges is imposed (by the 
pH of the solution). The positions of the charges along the chain are also not fixed; the 
charges can move by recombination and redissociation of an H + . The chemical potential 
of the charges n depends on the fraction / of dissociated acid groups and is related to the 
pH of the solution by [2] 

pH = pK + [x(f). (2.17) 

The charge chemical potential has two contributions, an entropic contribution related to 
the mixing of charged and non charged groups along the chains fcgTln[//(l — /)] and an 
electrostatic contribution jtt e /(/) = N~^dF e i/df which is the derivative of the electrostatic 
free energy of the chain. In a good or 9 solvent, the electrostatic contribution \i e \ is 
an increasing monotonic function of / and for a given pH the fraction of charges / is 
well defined, the properties of the chains are very similar to those of a quenched chain 
with a charge fraction /. In a poor solvent, the blob model introduced by Khokhlov and 
discussed above leads to a non monotonic variation of \i e \ and thus of /i with /. This non 
monotonic variation indicates a conformational transition of the chain between a collapsed 
weakly charged conformation and a stretched strongly charged conformation. The collapse 
transition of a polyelectrolyte with varying / is thus predicted to be a first order transition 
[25]. The discontinuous collapse of the chains could explain titration curves of polyacids 
or polybases in a poor solvent (polymethacrylicacid for example) which show a plateau in 
the pH — / curve [26] . 

3. Local aspects of screening 

The Flory-like model presented in the previous section considers only the interactions 
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between charged monomers along the polyelectrolyte chain and ignores the role of the 
small counterions that neutralize the polyelectrolyte charge. When the polyelectrolyte is 
strongly charged, the electrostatic potential on the chain is large and some of the counteri- 
ons remain bound to the chain. This phenomenon is known as counterion condensation, or 
Manning condensation [7]. For many purposes one must then consider that the chain has 
an effective charge due to the charged monomers and to the sheath of bound, or condensed, 
counterions. The effective charge is lower than the nominal charge of the monomers. Coun- 
terion condensation can be effective even at infinite dilution and reduces the electrostatic 
interactions between monomers. We first give here the more qualitative description of 
the condensation due to Manning [7] and then present an earlier discussion based on the 
Poisson-Boltzmann equation due to Fuoss et al. [27]. The condensed counterion sheath 
around a polyelectrolyte chain is electrically polarizable; this can induce attractive inter- 
actions between polyelectrolyte chains. We discuss in section 3.2. possible mechanisms for 
attractive electrostatic interactions in polyelectrolyte solutions. 

3.1 Counterion condensation 

For simplicity, we consider a negatively charged polyelectrolyte chain, in a salt-free 
solution, which has locally a rodlike conformation and can thus be considered as a rod 
over a distance L much larger than the monomer size b. The distance between charges 
along the chain is A and the fraction of charged monomers is / = b/A. If we ignore the 
contribution of the counterions, the electrostatic potential at a distance r « L from the 
chain is obtained from Gauss theorem as ip(r) = 2£ B /Aln(r) (we use here as unit of the 
electrostatic potential kgT/e where e is the elementary charge). The distribution of the 
counterions around the chain is given by Boltzmann statistics, for monovalent counterions, 
n(r) = no exp —ip(r) ~ r ~( 2 ^sM) where no is a constant. The total number of counterions 
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per unit length within a distance r from the chain is then 



r 



r 




2nr'dr'n(r') = j 




2 7 m r'( 1 - 2 V^V 



(3.1) 







When the charge parameter u = £ B /A is smaller than one, the integral giving p(r) is 
dominated by its upper bound and p(r) decreases to zero as the distance r gets small. In 
this case where the charge fraction is small f < b/£ B there is no counterion condensation 
on the polyelectrolyte chain. In the opposite limite where the charge is large u = £ B /A > 1, 
the integral giving p(r) diverges at its lower bound indicating a strong condensation of the 
counterions on the polyelectrolyte chains. As the condensation proceeds, the effective value 
of the charge parameter decreases. The counterion condensation stops when the effective 
charge parameter is equal to one. The polyelectrolyte chain and the condensed counterion 
sheath are then equivalent to a polyelectrolyte with a distance between charges A = £ B 
or to a fraction of charged monomers / = b/£ B . The remaining monomeric charges are 
neutralized by the condensed counterions. 

One of the directly measurable quantities that strongly depends on the condensation 
of the counterions is the osmotic pressure of the polyelectrolyte solution. As shown below, 
in many cases the osmotic pressure of a polyelectrolyte solution is dominated by the coun- 
terions, the contribution of the polyelectrolyte chains being only a very small correction. 
Below the Manning condensation threshold, the counterions are not condensed and are 
essentially free. At leading order, the solution can be considered as an ideal gas of coun- 
terions of concentration fc where c is the total monomer concentration and the osmotic 
pressure is [28] 



The first correction to the ideal gas behavior is due to the polarisation of the counterion gas 
by the polyelectrolyte chains. The polarisation energy is estimated using the Debye-Hiickel 
approximation for the gas of counterions [29] . The screening is due to the counterions and 



= k B Tfc 



(3.2) 
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k = 4n£ B fc. The polarisation free energy per unit volume reads 

F p oi = -k B Tfculn( K ) = -^ C (^|) ln(4n£ B cb/A) (3.3) 

The osmotic pressure of the solution is directly calculated from this relation 

tv = k B Tfc(l - u/2) (3.4) 

In many cases, for weakly charged polyelectrolytes, £ B ~ 6, so that u ~ /, and the 
polarisation correction is small. In most of the following, we use the ideal gas expression 
for the osmotic pressure of the counterions in a weakly charged polyelectrolyte solution. 

Above Manning condensation threshold u > 1, the bound counterions do not contribute 
to the osmotic pressure. The remaining counterions, at a concentration cb/l B = fc/u, 
behave as a Debye-Hiickel gas, polarized by a chain which has an effective charge parameter 
equal to unity. The osmotic pressure is given by equation (3.4), where fc is replaced by 
cb/£ B and u by 1, 

7T = k B Tcb/(2£ B ) = k B Tfc/(2u) . (3.5) 

Several other physical quantities such as the electrical conductivity of the solution or 
the electrophoretic mobility of the chains [30] strongly depend on the condensation of the 
counterions and are dominated by the free counterions. Their determination is, in general, 
in good agreement with Manning condensation theory. 

3.2 Poisson-Boltzmann approach 

Fuoss, Katchalsky, and Lifson [27] have studied the interaction between the polyelec- 
trolyte chains and their counterion clouds using the so-called Poisson-Boltzmann equation. 
They consider a solution of infinite rodlike molecules, which are all parallel. The number of 
molecules per unit area is Y. Around each of these molecules, there exists an equipotential 
surface where the electric field vanishes. On average, this surface can be approximated by a 
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cylinder parallel to the polyelectrolyte chains with a radius R defined such that TttR 2 = 1. 
The average density of counterions in the cylinder is q = fc = 1/(7tR 2 A) where A is the 
distance between charges along the chain. ^From an electrostatic point of view, each chain 
with its counterions in the cylinder of radius R is independent of all the other chains. The 
electrostatic problem that must be solved is therefore that of a single infinite chain with 
its counterions embedded in a cylinder of radius R with the boundary condition that the 
electric field vanishes on the surface of the cylinder. 

The electrostatic potential in this cell satisfies the Poisson equation Aip = —4n£ B Ci(r) 
where Cj(r) is the local counterion concentration at a distance r from the chain. The local 
counterion concentration is obtained from Boltzmann statistics Cj(r) = q exp (— ^(r)). 
This leads to the Poisson Boltzmann equation 

AV>(r) = -k 2 exp (-V>(r)) (3.6) 

The Debye screening length is defined here as k 2 = 4ir£ B cf, it is such that k 2 R 2 = 4u 
where u is the charge parameter defined above. The Poisson Boltzmann equation has been 
solved exactly in this geometry in reference [27]. The expression of the potential critically 
depends on the value of the charge parameter. Below the Manning condensation threshold 
where u < 1 the electrostatic potential is given by 

2 2 

^ = ln[^2- sinh 2 (-arctanh/3 + (3 ln(r/R)] (3.7) 

where (3 is an integration constant related to u by u = (1 — (3 2 )(1 — /3/tanh(/3 ln(6/i?)) _1 
Above the Manning condensation threshold the electrostatic potential reads 



ifi = In 



2 2 

sin (— arctan/3 + f3ln(r/R) 



(3.8) 



2f3 2 

and (3 is given by u = (1 — /3 2 )(1 — [3/ tan(/3 ln(6/i?)) _1 . The distribution of counterions 
around each polyion can then be determined using Boltzmann's statistics. 

18 



At each point of the solution, the pressure has two contributions, an ideal gas con- 
tribution and an electrostatic contribution proportional to the square of the local electric 
field. At the edge of the cell (r = R), the electric field vanishes and there is only an ideal 
gas contribution. The pressure being uniform throughout the solution can be calculated 
at this point it = k^Tc^r = R) = /c^Tqexp — (ip(r = R)). Both above the Manning 
condensation thereshold (3.4) and below the threshold (3.5) the calculated pressure is in 
agreement with that calculated from Manning theory. 

Manning theory describes the condensation of the counterions as a transition between 
two states, a bound state and a free state. In the Poisson Boltzmann approach, the dis- 
tribution of the counterions is continuous and there is no bound state. However, close to 
the polyelectrolyte chain, the interaction energy between one counterion and the polyelec- 
trolyte chain becomes much larger than kgT (in the region where ^ » 1). One could 
then divide space into two regions, a region close to the chain where the interaction energy 
is larger than k%T (with a size of the order of a few times the monomer size b) and where 
the counterions can be considered as bound and a region further away from the chain 
where the interaction between a counterion and the chain is smaller than kgT and where 
the counterions can be considered as unbound. 

3.3 Attractive electrostatic interactions 

The counterions condensed on a polyelectrolyte chain move essentially freely along 
the chain; they more or less form a one dimensional gas of average density g/b where 
g = f — b/£ B is the fraction of condensed counterions. The charge density along the chain 
is therefore not a frozen variable and shows thermal fluctuations due to the mobility of 
the counterions. When two chains are sufficiently close, the charge density fluctuations 
on the two chains are coupled by the electrostatic interactions; this leads to attractions 
between the polymers which are very similar in nature to the Van der Waals interactions 
between polarizable molecules. A fluctuation induced attraction can be expected for all 
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polyelectrolytes where the charges are mobile along the chain. This is for example the case 
above Manning condensation threshold but also for annealed polyelectrolytes, which are 
polyacids or polybases, where the charge can be monitored by tuning the pH of the solution 
[25? ]. The ionised groups are not fixed on these chains, only the chemical potential of the 
charges is imposed and the charge density also has thermal fluctuations. 

We first calculate the attractive interaction for two parallel rodlike molecules of length 
L at a distance i in a solution where the salt concentration is n, following the ideas 
introduced by Oosawa [31]. When a charge fluctuation 5c\(z) occurs along the first chain 
it polarizes the second chain where a charge fluctuation of opposite sign Sc2(z) appears. 
The attraction between the two chains is due to the Coulomb interaction between these 
charge fluctuations. Within a linear response approximation, the free energy of the charge 
fluctuations can be written in Fourier space as 

^^EE^M^H?) (3.9) 

9 i,j 

8ci(q) is the Fourier transform at a wave vector q of the charge fluctuation Sci(z) (i = 1,2). 
The charge structure factor on one of the chains is Sn(q) = g/b and corresponds to a one- 
dimensional ideal gas of charges. The crossed structure factor is due to the electrostatic 
interactions between the two chains obtained by summing the Debye-Hiickel interactions 

/ f ]2 u iqx cos 9 
4^2 fc 2 + g 2 + /c 2 = 2 ^0W<Z 2 + « 2 ) 1/2 ] (3-10) 

where Kq[x] is a modified Bessel function of the second kind [32]. The fluctuation con- 
tribution to the interaction free energy between the two rods is obtained by summing the 
partition function exp—H/kgT over all fluctuations. We obtain 

F a u(x) =L kj ff^^~ ^K Q[x ( q * + . 2 ) 1/2 ]) (3.11) 

This expression can be simplified by assuming that the fraction g of recondensed counte- 
rions is smaller than one and by looking for the asymptotic limits where the distance x 
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between the rods is smaller or larger than the screening length k . For short distances 
between the rods, the interaction energy is 

F a tt(x) = -*k B T 9 ^ (3-12) 
where a is a numerical constant. For large distances, the attractive interaction decays as 

F att (x) = -pk B T 9 ^ £ exp(-2^) (3.13) 

(3 being another numerical constant. As expected the interaction is attractive. It is propor- 
tional to the square of the Bjerrum length indicating that this is a second order electrostatic 
effect: the fluctuation on the first rod creates an electric field that polarises the second rod; 
the charge fluctuation on the second rod then creates an electric field which interacts with 
the charge fluctuation of the first rod. For parallel rods, the interaction is also proportional 
to the length of the rod, each segment interacts mainly with the segment of the other rod 
which is directly facing it. For the same reason, the attractive interaction is screened over 
a length k _1 /2, i.e. over half the screening length of the repulsive coulombic interaction. 

When the two rods make an angle 6 (0 < 6 < tt/2), the attractive electrostatic inter- 
action cannot be calculated from the fluctuation free energy (3.10). Instead, we directly 
discuss it in terms of the polarisability of the two polyelectrolytes 2 If a charge fluctua- 
tion Sci(q) is created on the first rod, it creates an electrostatic potential on the second 
rod 5ip(q') = K(q, q')5c\(q) where the kernel is defined as K(q, q') = exp ^~ px ^ with a 
wavevector p related to q and q' by p 2 = k 2 + q +q gin 2 g g C ° S ^ ■ The response function of the 
one dimensional gas of recondensed counterions on the second rod is | and the induced 
charge fluctuation on the second rod is 8c2(q') = —^K(q, q')Sci(q) . This charge fluctua- 
tion creates a potential on the first rod S^'(q) = K(q, q')Sc2(q')- The attractive interaction 

This is equivalent to performing a perturbation expansion to second order in the crossed 
term that appears in equation (3.10). 
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energy between the two rods forming an angle 9 at a distance x is then estimated as 

F att (x) = -k B T 9 - J d -^-K\q, q>) < 8c\{ q ) > (3.14) 

The average value of the charge fluctuation is given by the ideal gas statistics < 5cf(q) > = 
g/b. This leads to an attractive interaction between non parallel rods 

2 2 

F att (x) = -k B T-^-(n - 0)E i (2Kx) (3-15) 
o z sin 9 

where Ei(u) is the exponential integral function [32]. As for parallel rods, the attractive 
interaction is a second order electrostatic effect, it is proportional to l B and it decays at 
large distances as exp —2kx. When the angle between the two rods is not very small, 
the attractive force is independent of the length of the rods. In this geometry only the 
chain segments in the vicinity of the crossing point interact. If the distance between the 
rods is smaller than the screening length, the crossover to the parallel rod behavior occurs 
when the angle 9 is smaller than x/L. At short distances x < k -1 , Ei(tt) = ln-u and the 
attractive interaction energy varies logarithmically. As the attraction is independent of 
the chain length, the attractive interaction between two polymers which are rodlike over 
a distance larger than the screening length is also given by equation (3.14). 

The fluctuation attractive interaction must be compared to the repulsive coulombic 
interaction between the chains. For two rods at a distance x that make a finite angle, 
the repulsive Coulomb interaction is of order Fel ~ k B T2n£ B f' 2 '/(Kb 2 ')exp—Kx where 
/ is the effective fraction of charged monomers (equal to bji B above the condensation 
threshold). At large distances, as already mentionned, the screening of the attractive 
interaction is stronger and the repulsive interaction is dominant. At small distances the 
attractive interaction is dominant. The distance at which the attractive force is larger 
than the repulsive force is of order x ~ £ B (g/ f) 2 . In general this corresponds to very short 
distances of the order of the Bjerrum length. The attractive fluctuation interaction is thus 
expected to be important only in strong coupling situations where the Bjerrum length 
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is large (for example in the presence of multivalent ions) or if some external constraint 
imposes very small distances between chains. 

The attractive interaction can also be considered within the framework of the Poisson- 
Boltzmann approach. The Poisson-Boltzmann equation is derived from a mean field theory 
that considers only the average concentration profile of the counterions around each poly- 
electrolyte chain. A more refined theory should take into account the thermal fluctuations 
around this concentration profile. Such a theory has been built to study the interactions 
between charged colloidal particles and also leads to attractive interactions [33,34]. Quali- 
tatively, the conclusions are the same as the one presented here, the attractive interactions 
are dominant at short distances and becomes relevant only in the limit of strong electro- 
static coupling. 

Experimentally, the addition of multivalent ions to polyelectrolyte solutions often pro- 
vokes a precipitation. This can be explained by attractive electrostatic interactions [35] or 
by complexation of the polymer by the multivalent ions [36] . 

The interaction between charge fluctuations is not the only mechanism that leads 
to attractive interactions between polyelectrolyte chains. Recently, Ray and Manning 
[37] have presented a model for attractive interactions based on the overlap between the 
condensed counterions sheaths around the two interacting polymers. In the simple version 
of Manning condensation theory presented here, the polyelectrolyte chains are infinitely 
thin lines and the pointlike counterions condense directly on the line. A more refined 
version allows for a finite condensation volume that can be determined by minimization 
of the free energy of the condensed counterions. When two parallel rods come close to 
one another, the condensation volume around each rod expands in the space between 
the two rods and the two condensation volumes overlap. This leads to the formation of a 
polyelectrolyte dimer with a single condensation sheath. The expansion of the condensation 
volume leads to an increase in translational entropy of the counterions and thus to an 
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attractive interaction. The model shares thus some analogy with the capillary condensation 
between thin films wetting parallel cylinders [38]. An attractive interaction is predicted 
at intermediate distances much larger than the polyelectrolyte diameter but much smaller 
than the electrostatic screening length. This interaction has been calculated only for 
parallel polyelectrolyte rods, it has not been calculated for rods intersecting at a finite 
angle. 

4. Electrostatic rigidity 

In this section, we assume that the interaction between charged monomers can be 
described by the Debye-Hiickel potential, (2.16). Within this approximation, we discuss 
the structure of the polymer chain at length scales intermediate between the short length 
scales discussed in section 3. and the chain size. It has already been mentioned in section 
2.4 that at very large scales, a charged chain has the same structure as a neutral polymer 
chain in a good solvent, since the potential (2.16) is short ranged. The radius increases with 
molecular weight as R ~ iV 3 / 5 (in the Flory approximation). At short scales (L < k _1 ), 
on the other hand, screening is inoperant, so that the chain is expected to take the rodlike 
structure described in section 1. The simplest assumption concerning the chain structure 
at intermediate length scales, which was made in all early work in the field [39] , is that the 
crossover from rodlike to Flory-like behaviour takes place above a length scale k -1 , equal 
to the range of the interactions. This simple assumption was challenged by the works of 
Odijk [8,10] and Skolnick and Fixman [9], who showed how the Debye-Hiickel interaction 
can induce a rodlike conformation at length scales much larger that the interaction range. 
Their theory, and its limitations, are presented here. 

4-1 The Odijk- Skolnick- Fixman theory 

In their calculations, Odijk-Skolnick and Fixman (OSF) consider a semi-flexible chain 
with total contour length L characterized by its "bare" persistence length £ Q , and that 
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carries charges separated by a distance A along its contour. The interaction between the 
charged monomers is given by (2.16), and only electrostatic interactions are taken into 
account (the polymer backbone is in a theta solvent). For strongly charged chains with 
A < £ B Manning condensation can in a first approximation be accounted for by replacing 
A with £ B . In the limit £ << A of a weakly charged chain, this model is equivalent to the 
gaussian chain model discussed in section (2.1), with N = L/A, b 2 = A£ Q and / = 1. We 
consider here the general case, where £ /A is not small [40]. 

The total energy of the chain is the sum of the intrinsic curvature energy and of the 
(screened) electrostatic energy. If the polymer is described as a planar curve of curvilinear 
length L, a chain configuration is specified by the function 9(s) defined for —L/2 < s < L/2 
by cos(0(s)) = t(s) • t(0) where t(s) is the unit vector tangent to the chain at the point 
of curvilinear abscissa s. Assuming that 9(s) remains small, the electrostatic potential 
between two charges located at s± and s 2 on the chain may be expanded around the value 
obtained for a rodlike configuration, vdh(\s2 — si |) - The following expression is then 
obtained for the total energy H[9] of a given configuration {9(s)}: 

H[9] =H + \k B T [ L/2 ds [ L/2 ds' ^ [t S(s - «') + K{s, s')} ^ . (4.1) 

2 J-L/2 J-L/2 ds ds 

Here Hq is the electrostatic energy of a rod, the term proportional to £ is the "bare" curva- 
ture energy of a noninteracting semiflexible chain, and the contribution of the electrostatic 
interactions is described by the kernel K(s, s') which reads (for s > s'): 

K(s,s') = ^[ S ds, [ L/2 ds 2 v' DH (s 2 -s 1 ) {S2 - s){s/ - Sl) (4.2) 
A J-L/2 Js {$2 - 8i) 

with v' DH (s) = dv ^ s H The only assumption made in obtaining this result is that, within the 
range of the interaction potential, the chain remains in an almost rodlike configuration, 
i.e. \s2 — si| — — R| << |s2 — s l| f° r \ s 2~ s l\ < K -1 - For long chains, the integration 

in (4.2) can be extended to infinity, and K(s, s') becomes a function of s — s', 

i roo 3 
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The statistical properties of the chain are obtained by integrating the Boltzmann factor 
exp(— H[6]/kBT) over all possible configurations, i.e. over all functions 0(s). The inte- 
gration can be carried out analytically, since the energy is a quadratic function of 9. The 
calculation is simplified by introducing the Fourier transform of the kernel (4.3), 

rOO 

K(q) = / dsexp(iqs)K(s) 



^OSF 9 1 2 M 9 

q z \ Q V K 

where 

W=^/(4A 2 « 2 ) (4.5) 

has the dimension of a length, and is known as the Odijk-Skolnick-Fixman length. The 
local flexibility of the chain can be characterized by the mean squared angle < 6>(s) 2 > 
between the chain direction at the origin and after a contour length s. For a neutral 
semiflexible chain, < #(s) 2 >= s/£ varies linearly with s. The chain configuration can be 
described as resulting from an "angular diffusion" process, with a diffusion constant £~^. 
For the charged chain, with the approximate energy (4.1), one gets 



< 6(sf >=- I dq — "J 1 "' V— , (4.6) 



4 r°° sin 2 (gg/2) 1_ 

7T JO Q Q 1 £ Q + K(q) 



This expression simplifies in the limits of small s and large s. For large s, 

< 0(s) 2 >= T ^ J — (4.7) 

which expresses the fact that at large scales, the chain conformation can be described by an 
effective persistence length £ + £ OSF: which is the sum of a "bare" and of an electrostatic 
contribution. This is the well known result first obtained in references [8] and [9]. It 
indicates that the influence of the screened electrostatic interactions can extend much 
beyond their range since l OSF is, for weakly screened solutions, much larger than 

The persistence length also decreases with the salt concentration as t OSF — n^ 1 
whereas the Debye screening length has a slower decay k^ 1 ~ n -1 / 2 . 
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For small values of s, (4.6) reduces to 

< d(sf >= (4.8) 

The chain statistics at short scales are not modified by electrostatic interactions. The 
crossover between the "intrinsic" regime described by (4.8) and the "electrostatic" regime 
described by (4.7) takes place when the electrostatic interactions become strong enough to 
perturb the statistics of the neutral flexible chain. The crossover length s c can be obtained 
qualitatively from the following argument If a small chain section, of length s < ft -1 , is 
bent to form an angle 6>, the cost in "bare" curvature energy is fcgT£ o 2 /s, while the cost 
in electrostatic energy is k B T£ B (s/A) 2 (9 2 / s). The two energies are comparable for s ~ s c , 
which gives s c ~ A(£ /£ B ) 1 / 2 . 

A more detailed treatment of the crossover regime is possible by ignoring logarith- 
mic factors and replacing the exact Kernel (4.4) by the approximate expression, K(q) = 
£ osf k 2 /(k 2 + 2q 2 ). The integral (4.6) can then be computed analytically, and yields, 



< 9(8)2 >= 1 ~TT + 1/2 — TtJ 1 ~ ex P("f )) (4.9) 

K £ 1 J 2 (£ o+ £ osF )3/'2 Sc 



with s c = k 1 (£ OSF +£ )' This formula crosses over from the "intrinsic" regime 

towards the electrostatic regime for s ~ s c . In the weak screening limit, £ << £ OSF and 
the simple result s c ~ A(£ /£ B ) 1 / 2 is recovered. 

The picture that emerges from this calculation is that the chain flexibility depends on 
the length scale. At short scales, s < s c , the chain structure is determined by its bare 
rigidity £ , while the electrostatic rigidity (4.5) dominates at large scales. Large scale 
properties, such as the giration radius or the structure factor at small wavevector, can be 
determined by applying the standard formula for semi-flexible chains of persistence length 
£ OSF - Excluded volume effects between Kuhn segments of length 2£ OSF are accounted 



3 From its definition, s c can be described as the contour length of an "orientational blob" 
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for by assigning a diameter k to each segment. The excluded volume between Kuhn 
segments is £^ sf k~ 1 [10]. 

The only approximation required to obtain the Odijk-Skolnick-Fixman length is the 
expansion that yields equation (4.1). The calculation is therefore consistent if the angle 
< 6(k~ 1 )^ > is small compared to unity. It is easily checked that for s ~ the second 
term dominates in the r.h.s. of (4.9). The requirement < 6>(k -1 ) 2 ><< 1 is then equivalent 
to s c £ OSF /(£ Q (£ Q + £ OSF )) — s c /£ Q << 1. i-e. the angular deflection of the chain must be 
small when the crossover region is reached. In other words, the angular fluctuations that 
take place before the electrostatic interactions can come into play and rigidify the chain 
should not be too strong. If these fluctuations are too large, i.e. if the chain is too flexible, 
the perturbation expansion that underlies the OSF calculation breaks down. The criterion 
for the validity of the calculation can be simply written, in the limit of weak screening, 
as £ Q > A^/£ B . This implies that the OSF calculation should be directly applicable to 
stiff chains such as DNA (£ ~ 50nm), but has to be reconsidered for flexible chains such 
as polystyrene sulfonate (£ ~ lnm). Indeed, numerical simulations clearly confirm the 
behaviour described by equations (4.4) and (4.6) for values of s c /£ Q smaller than 0.2 [40]. 
Deviations from this behaviour appear for larger values of s c /£ . 

Before closing this section, it is worth mentioning that the electrostatic rigidity of a 
charged chain was computed numerically by Le Bret and Fiman [41,42] using the Poisson- 
Boltzmann equation, rather than the Debye-Hiickel approximation. At high ionic strength, 
the results deviate significantly from the OSF prediction, and the electrostatic rigidity 
tends to behave as rather than 

4-2 Alternative calculations for flexible chains 

The calculation presented in the previous section shows that the Odijk-Skolnick-Fixman 
theory in its original form is not applicable to flexible or weakly charged chains as such. 
A generalization of this theory was nevertheless used by Khokhlov and Katchaturian [43] . 
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These authors propose that the expression (4.5) can be used for flexible chains, with the 
only modification that the distance between charges A had to be replaced by £ e /(/<?) (see 
equation (2.7). The bare persistence length £ is also replaced by the blob size £ e (A recent 
calculation of Li and Witten, including the thermal fluctuations of the chain gives a more 
quantitative basis to this approach [44]). As a result, the total persistence length of the 
chain reads 

e '"^^ + ik- (4 - 10) 

The arguments of section (4.3) show, however, that the extension of the OSF formula is 
rather speculative, since the condition that was derived for the validity of the OSF theory 
is not met for a chain characterized by the parameters £ = £ e and A = £ e / fg. When the 
persistence length is dominated by the electrostatic interactions, the radius of the chain is 
given by the Flory statistics 

RKK = N^i-^K-^b^ . (4-11) 

This generalisation of the Odijk theory is consistent if the Debye screening length is larger 
than the electrostatic blob size £ e - At higher ionic strength, in the weak coupling limit 
there is no electrostatic rigidity and the polyelectrolyte chain behaves as a flexible chain 
with short range interactions given by the Debye-Hiickel formula. Its radius is given by 
equation (4.12) below. 

A number of approximate theories, based on different variational ansatz, have been 
proposed to describe the structure of flexible chains [45,40,46,47]. These calculations differ 
at the technical level. The basic idea is to describe a flexible charged chain (with screened 
interactions) by some model of noninteracting semiflexible chain, and to variationally opti- 
mize the persistence length of the noninteracting system. All these variational calculations 
produce for flexible chains a persistence length that scales as i.e. the "naive" result 
of Katchalsky and Pfeuty. An exemple of such a variational calculation can be given for 
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the particularly simple case of a freely jointed charged chain, with £ = A. The variational 
ansatz used to represent the system is that of a noninteracting chain in which the angle 
between neighboring segments is bounded by a maximum value 9 ma x, so that the persis- 
tence length is £ p = 2A/9^ ax . The variational entropy per segment is then — ln(9 max ). 
The variational energy can be expressed as a function of the structure factor of a chain 
with a persistence length £ p , S n i(q) as J d^c[S n i(q)vjjjj(q). Using the expression of the 
structure factor obtained by desCloizeaux [48], the energy is the sum of the electrostatic 
energy of a rod and of a correction due to the bending that can be approximated by 
kBT{£ B /A){l/K£p). Minimizing the sum of these two terms with respect to 6 max , one 
finds £ p ~ This approach can be straightforwardly generalized to the case where each 
rod in the freely jointed chain represents an electrostatic blob [40] , in which case it applies 
to the charged gaussian chain considered in section 2. The variational model is then that 
of a of a " semiflexible chain of electrostatic blobs" , analogous to the " chain under tension" 
model (section 2), except for the fact that the direction of the tension is now fluctuating, 
and becomes uncorrelated over a persistence length k~ 1 . In that case, it must be realized 
that the "persistence length" is defined in reference to the distance in space, rather 
than the contour length along the chain. The radius of the chain scales in this case as 

R P = iV 3 / 5 £- 3 / V 2 / 5 6 6 / 5 . (4-12) 

Variational calculations have the obvious drawback that the result crucially depends 
on the variational ansatz. An ansatz too far from the actual structure can yield incorrect 
results. Also, their range of validity is difficult to assess. The variational calculation for 
gaussian chains briefly described above is consistent only as long as the "electrostatic blob" 
concept can be used, This gives an upper limit for the charge density of the chain for which 
its description by a "semiflexible chain of electrostatic blobs" is possible, A -1 < (£ £ B ) ~ 1//2 
(this corresponds to g > 1 in 2.7). This upper limit corresponds to the lower charge 
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density for which the OSF calculation is expected to be consistent. The coincidence, 
however, might be only fortunate. A complete calculation that would crossover from the 
behaviour of the electrostatic stiffness which is well understood for rigid chains, to the 
k^ 1 behaviour that is suspected for flexible chains on the basis of variational calculations, 
is still missing. 

4-3 The case of poor solvents 

The above discussion focuses on the case where the neutral chain backbone is in a 
good or theta solvent. A rather different result is obtained for the case of a flexible, weakly 
charged chain in a poor solvent. This case was considered in the absence of screening in 
section (2.1), where the chain was described as a linear string of "poor solvent blobs" , each 
of size £ e = Hf 2 ^/^) 1 ^ an d containing rf(f 2 £ B /b)~ 1 charged monomers (cf. equation 
(2.10). Following the procedure of [43], one can attempt to apply the OSF theory to 
the chain of blobs. The theory is valid if the linear charge density along the chain of 
blobs, fg/ie, is larger than the critical value (ie£ B )~ 1 ^- This condition can be rewritten 
as r > (/ 2 £ B /6) 1 / 3 . For weakly charged chains (/ << 1), it is thus satisfied even for 
moderately poor solvents. The electrostatic persistence length that results from applying 
the OSF theory is 



The surprising implication of this result is that the persistence length, (and therefore the 
radius of gyration of the chain), increases as the solvent quality decreases. This increase 
is the consequence of an increase in the linear charge density, which is itself induced 
by the collapse of the monomers into dense blobs. This effect, however, exists only as 



increase in the persistence length as r increases should occur in the temperature range 



— £,e + 




(4.13) 



long as the charge density exceeds the value t B at which Manning condensation takes 
place. This gives an upper limit for r, r < f(f^£ B /b)~ l l^. The conclusion is that an 
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5. Charged Gels and Brushes 

Many properties of polyelectrolyte solutions are dominated not by the chain confor- 
mation, but by the counterions. We discuss below two examples of this situation, the 
polyelectrolyte gel and the polyelectrolyte brush. In both problems, the polyelectrolyte 
chains occupy a small region of space (close to the grafting surface for the brush), sur- 
rounded by pure solvent. The small ions in the solution are free to diffuse in and out the 
polyelectrolyte region, and the external solvent region acts as a reservoir that imposes the 
small ion chemical potential. A Donnan equilibrium [49,50] is reached, that relates the 
concentrations of the small ions inside the polyelectrolyte region to their imposed concen- 
tration in the reservoir. If the monomer concentration in the polyelectrolyte region is c, 
the charge density due to the polymer is fc. The reservoir contains salt at a density n, and 
the Debye-Hiickel screening length is given by = 8nn£ B . For simplicity we assume 
in the following that the counterion of the polymer is identical to the positive ion of the 
salt, we call n+ and the concentration of positive and negative ions respectively in the 
polyelectrolyte region. 

The polyelectrolyte charge creates an electrostatic potential difference between the 
polyelectrolyte region and the reservoir U. The chemical potentials of the small ions in the 
polyelectrolyte region are /i+ = /c^Tlnn-i- + U and \i- = k^Tlnn^ — U. At equilibrium 
these chemical potentials are equal to the salt chemical potential in the reservoir /j = 
kgTlnn. This leads to a chemical equilibrium law for the small ions 

n+n- = r? (5-1) 

A second relation between the concentrations n+ and n_ is provided by the electroneu- 
trality condition n+ = n_ + fc. The salt concentration inside the polyelectrolyte region 
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(given by is smaller than that of the reservoir. The difference in osmotic pressure 
between the two regions is n = k B T(n+ + n-- 2n) = [(/c) 2 + An 2 } 1 / 2 - 2n. When the 
salt concentration n is smaller than the counterion concentration fc, the difference in salt 
concentration is small and the osmotic pressure is that of the ideal gas of counterions 
discussed in section (3.1), tt = kgTfc. In the limit where the salt concentration is larger 
than the counterion concentration, the osmotic pressure is given by 

This osmotic pressure defines an effective virial coefficient between monomers v e \ = f 2 /2n = 
4ir£ B f 2 K~ 2 . This is identical to the excluded volume that can be calculated from the 
Debye-Hiickel interaction (2.16) between monomers. 

The major assumption made here is that the ion densities are uniform. This is the case 
when regions of size k" 1 around each polyelectrolyte, where the counterion concentration 
is increased, overlap, i.e. when the Debye-Hiickel screening length is larger than the dis- 
tance between chains. This is always true in the absence of salt. The osmotic pressure 
calculated from the Donnan equilibrium is of purely entropic origin. If the ion densities 
are uniform the electrostatic contribution to the osmotic pressure is small compared to the 
translational entropy of the small ions as checked below for a gel. The osmotic pressure 
is thus independent of the strength of the electrostatic interaction characterised by the 
Bjerrum length £ B . The only role of the electrostatic interaction is to enforce electrical 
neutrality. 

5. 1 Grafted polyelectrolyte layers 

Grafted polymer layers (polymer brushes) [51] are obtained by attaching one of the 
polymer end points on a planar solid surface. If the distance between chains is smaller 
than their natural size, a thick layer with a size proportional to the chain molecular weight 
is formed. The chains are stretched by the repulsive interactions between monomers. 
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Polymer brushes are often very efficient to enhance colloidal stabilisation [52] . Grafted 
polyelectrolyte layers have been described theoretically by several authors [53,54,55,56] 
we follow here the lines of the original work of Pincus which assumes that the polymer 
concentration is uniform throughout the thickness h of the brush. A similar approximation 
for neutral polymer brushes was introduced by Alexander [57] and de Gennes [58]. We 
also assume in this section that the local electrostatic screening length is larger than the 
distance between chains imposed by the distance D between grafting points so that the 
counterion concentration is roughly constant in the layer. 

The electrostatic field and the distribution of counterions in the vicinity of a charged 
solid surface with a surface charge density p can be determined from the Poisson-Boltzmann 
equation. In the absence of added salt, the counterion density decays as a power law of the 
distance from the surface, and the counterions are confined in the vicinity of the surface 
over the Gouy-Chapman length A = (2np£ B )^ 1 . If the grafting surface is neutral, the 
charge density of a grafted polyelectrolyte with a = D~ 2 grafted chains per unit area layer 
is due to the monomer charges p = crNf, the Gouy-Chapman length [23] of the polymer 
brush is 

A * ^Tm B (5 ' 3) 

Two cases must then be considered. When the Gouy-Chapman length is larger than the 
brush size h, most of the counterions are outside the brush, and the brush is charged. 
When A is smaller than the brush height, the counterions are confined in the grafted layer 
and the brush is neutral. 

In the neutral limit, the brush is swollen by the osmotic pressure of the counterions 
tv = kgTfc. This pressure is balanced by the elasticity of the polymer chains. For a 
gaussian chain, the elastic pressure is 7r e / = ah/(Nb 2 ). The monomer concentration in the 
brush being c = a N/ h the thickness of the brush is given by 

h = Nbf 1 ' 2 (5.4) 
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The thickness of the brush is in this regime independent of the strength of the electrostatic 
intractions since the stretching of the chains is due only to entropic effects. The thickness 
given by equation (5.4) is larger then the size of an individual chain in a dilute solution 
given by (2.5). This result has been obtained by assuming that the monomer concentration 
is constant inside the adsorbed polymer layer. Similar scaling results are obtained if this 
constaint is released. Using a self-consistent field method, Zhulina et al. [54] have shown 
that the concentration profile has then a Gaussian decay. 

In the charged limit, one can in a first approximation assume that the counter ion 
concentration is constant up to the Gouy-Chapman length, it is equal to q = (fch)/X. 
The counterion pressure that stretches the chain is therefore 7r = kgTfch/X. The thickness 
of the brush is then 



In this regime the electrostatic interactions have a contribution of the same order as the 
entropic contribution and the thickness depends on their strength. The thickness of the 
layer grows faster than linearly with molecular weight, however one should keep in mind 
that the layer thickness remains smaller than the Gouy-Chapman length and thus that 
this thickness remains smaller than the thickness in the neutral regime (5.4). A precise 
determination of the counterion profile in this regime has been done by Pincus [53], the 
same scaling laws are obtained when this profile is taken into account. 

When salt is added to the layer, the polyelectrolyte brush behaves as a neutral polymer 
brush with an effective excluded volume v e \ = / 2 /(2n); the thickness of the brush is then 



The monomer concentration profile has also been calculated in this limit by Zhulina et al. 
[54] The important result is that the cross-over to the neutral brush regime (5.4) occurs 
when the salt concentration is larger than the counterion concentration in the brush. This 



h = N 3 £ B b 2 af 2 



(5.5) 




(5.6) 
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leads in general to a high ionic strength and the thickness of a grafted polymer layer is 
rather insensitive to ionic strength over a broad range of ionic strength. This makes grafted 
polyelectrolyte layers particularly interesting to promote colloidal stabilisation. 

If the polyelectrolyte is annealed, the pH inside the grafted layer can be significantly 
different from outside. This could induce a non monotonic variation of the thickness with 
grafting density [59]. 

5.2 Polyelectrolyte gels 

A model based on the Donnan equilibrium, very similar to that of polymer brushes, 
can be made for charged polymeric gels [60]. The polymer gel is at equilibrium with a 
reservoir of solvent if the osmotic pressure in the gel is equal to the osmotic pressure of 
the salt in the reservoir, i.e. if the swelling osmotic pressure due to the entropy of the 
counterions is balanced by the pressure due to the elasticity of the polymer chains. If the 
chains between the crosslinks of the gel are gaussian chains and have N monomers and an 
end to end distance (mesh size) R, the elastic pressure is of order n ~ k B T(c/N)(R 2 /Nb 2 ). 
If we now make the c* assumption proposed by de Gennes [4] that, in the swollen gel at 
equilibrium, the chains are just at the overlap concentration, the mesh size of the gel is 
such that c = N/R^. In the absence of salt, the mesh size is then given by equation (5.4) 

R = Nbf 1 ' 2 (5.7) 

In the presence of salt, the mesh size is given by 

i ? = iV 3/5 6( Z! ) l/5 (5.8) 
2n 

The important point to notice is that, in contrast to neutral gels, the mesh size of the 
gel is different from the radius of an isolated polymer chain in a dilute solution. As 
explained above, these results are valid when the salt concentration is small enough that 
the Debye-Hiickel screening length is larger than the gel mesh size i.e. when the mesh size 
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is larger than the radius of an isolated chain in solution given by (2.4). The chains in this 
case should be viewed as chains under strong tension. The small ions osmotic pressure is 
exerted at the surface of the gel and the force is transmitted to the internal chains via the 
crosslinks. 

The mesh size and thus the equilibrium swelling of the gel do not depend on the 
strength of the electrostatic interaction £ B . When the monomer concentration in the gel is 
uniform, the electrostatic energy of the gel vanishes. The actual concentration in the gel 
is not homogeneous and the inhomogeneity can be characterised by the structure factor 
S(q). The electrostatic energy per unit volume of the gel can then be expressed as 

E e i = c J - / ddS(q)^^ (5.9) 

The structure factor of the charged gel is not known exactly but one can expect it to have 

a sharp peak at a position q* ~ R^ 1 the height of the peak being of order N. The integral 

giving the electrostatic energy of the gel is dominated by the peak of the structure factor 

and can be estimated up to numerical prefactors 

ck B T fN 2 £ B 1 , , 

Ed = ~N R—T^tf (5 - 10) 

It can then be direcly checked that, in the weak screening limit where kR « 1, the 
electrostatic energy is smaller than the translational entropy of the counterions. 

In the strong screening limit kR » 1, the electrostatic interactions are dominant. The 
charged monomers in the gel interact via a Debye-Hiickel potential and the electrostatic 
interactions contribute to the rigidity of the chains between crosslinks as discussed in the 

previous section. At very high ionic strength the electrostatic persistence length is small, 

f 2 

and the chains behave as neutral chains with an electrostatic excluded volume v e \ = fa. 
The same scaling result as equation (5.8) is then found for the mesh size. 

The elastic shear modulus of the gel can be calculated by imposing a uniform defor- 
mation 7 to the crosslinks. The increase of the square of the mesh size imposed by the 
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deformation is of order (7-R) 2 and if the chains are gaussian, the elastic energy stored in 
the gel is 

F, = kr,T-^ (5.11) 
el KB N Nb 2 J 

This leads to a shear modulus 

c R 2 

G = h^T— — -tt ~ kgTfc (salt free case) 

NNb2 v (5.11) 

~ k B T (with added salt) 

n 

In all cases, the shear modulus of the gel is thus of the order of the osmotic pressure of 
the small ions (counterions and salt) 

The last property of the gel that we discuss is the cooperative diffusion constant, that 
can be obtained from a two-fluid model [61]. The displacement field of the gel w(r, t) results 
from a balance between the elastic restoring force and a viscous force due to the solvent 
drag. The elastic force per unit volume of the gel is GV 2 «. The viscous force exerted by a 
solvent flowing through a porous medium of pore size R is given by the Brinkman equation 
[62] and is of order r]R~ ^ where r\ is the solvent viscosity. The force balance on the gel 
is then written as 

V R~ 2 ^ = GV 2 u (5.12) 
ot 

The relaxation of the gel deformation is diffusive with a cooperative diffusion constant 

D = — (5.13) 
V 

The cooperative diffusion constant is independent of the molecular weight in the absence 
of salt (D ~ f 1 / 2 ). It decreases with iV in the presence of salt (D ~ A?" _2 / 5 (/ 2 /n) 1 / 5 ). 

Experimentally, the properties of charged gels, as well as the properties of neutral gels, 
strongly depend on the preparation condition of the gels. The simple results presented 
here (based on the so-called c* theorem) implicitly assume that the gel was prepared by 
crosslinking chains in the vicinity of their overlap concentration. A more refined theory 
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[63] based on an affine deformation hypothesis first proposed by Flory [64] and then de- 
velopped by Panyukov [65] has recently been constructed. It provides an explanation for 
the experiments of Candau and coworkers [66] which show that the elastic shear modulus 
of polyelectrolyte gels with a fixed concentration smaller than the equilibrium swelling 
concentration decreases with the fraction of charged monomers. 

6. Semidilute solutions 

The overlap concentration of polyelectrolyte solutions in the absence of salt is very low; 
the chain radius given by equation (2.4) increases linearly with molecular weight and the 
overlap concentration decreases as c* ~ b 4: /(£ B N 2 f' 2 ). Most of the experiments with long 
chains are thus performed in the semidilute range where the electrostatic interactions be- 
tween different chains are strong. Our understanding of the conformation of polyelectrolyte 
chains in semidilute solutions is however rather poor and no general view is available. We 
discuss in this section a few aspects of the static properties of interacting polyelectrolyte 
chains. 

6. 1 Ordering transitions in polyelectrolyte solutions 

The models presented in section 2 to describe the conformation of polyelectrolytes in 
salt free dilute solutions consider only a single chain and ignore the screening due to the 
counterions. If the polymer concentration is finite, the counterion concentration is finite 
and the screening length due to the counterions is given by = 4tc£ B fc. In a dilute 
solution the screening length is always larger than the average distance between chains 
d ~ (iV/c) 1 /^. Different chains thus interact via a long range pure Coulomb potential. 
The polyelectrolyte chains are strongly charged objects and this long range interaction can 
lead to the formation of Wigner crystals, very similar to the colloidal crystals observed in 
solutions of charged spherical particles [67] or charged elongated particles (viruses) with a 
mesoscopic size. The centers of mass of the particles are then regularly distributed on a 
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periodic lattice. The solution has a tendency to crystallise when the interaction between 
neighboring chains kgTN 2 f 2 £ B /d is larger than the thermal excitation kgT. This naive 
estimate predicts crystal formation at a very low concentration, way in the dilute range 
c ~ l/(N^f®£ B ). The crystal is expected to melt when screening is important, in the 
vicinity of the overlap concentration. There is no clear experimental evidence for this 
crystallisation in flexible polyelectrolyte solutions. This may be due to the weak elastic 
resistance of this crystal that could be destroyed by any kind of perturbation (mechanical 
perturbation) . In any case we expect strong interactions in dilute polyelectrolyte solutions 
that lead to a structuration of the solution and thus to the appearance of a peak in the 
structure factor S(q) of the solution corresponding either to liquidlike or to solidlike order. 
The wavevector q* at the peak is of the order of the inverse distance between chains 
q* ~ 27i/d ~ (c/iV) 1/3 [68]. The peak is observed experimentally and the variation of its 
position with concentration is in good agreement with this prediction. It can be suppressed 
by adding salt to the solution, so that the screening length k~ 1 becomes smaller than the 
distance between chains. 

As in dilute solutions, polyelectrolytes have, in a semidilute solution, a locally rodlike 
conformation that can be characterized by a persistence length £ p . Quite similarly to rod- 
like molecules, semiflexible macromolecules undergo an Onsager transition [69,70] between 
an isotropic liquid phase and a nematic ordered phase where the molecules are parallel. The 
Onsager concentration where this transition takes place has been calculated by Semenov 
and Khokhlov [71,72], for neutral molecules of diameter d it scales as c ~ l/(£ p db). For 
charged molecules, the diameter of the molecule must be replaced by an effective diameter 
equal to the range of the electrostatic interaction k~ 1 [70] and the Onsager concentration 
is of order c ~ l/(£ p K~ l b). All the models for polyelectrolyte conformation based on 
the Odijk-Skolnick-Fixman approach predict an isotropic-nematic transition for polyelec- 
trolyte solutions in the presence of salt. Experimentally, the transition is only observed 
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for very rigid polymers (such as DNA where a cholesteric phase is observed) where the 
intrinsic persistence length is very large and dominates the electrostatic persistence length. 
There is no definite evidence for an Onsager transition in flexible polyelectrolyte solutions. 
This is an unresolved issue. Possible explanations are based on the idea that the scaling 
theories ignore some prefactors that may be large and shift the transition to an unobserv- 
able value. The isotropic-nematic transition is also expected only for anisotropic enough 
objects, the ratio between the persistence length and effective diameter k" 1 must be larger 
than a finite number ~ 5 [73] . This criterion is not always met experimentally. A thorough 
theoretical study of possible isotropic-nematic transition in polyelectrolyte solutions has 
been made by Nyrkova [74]. In the following we will ignore this possibility and consider 
that a semidilute polyelectrolyte solution remains isotropic at any concentration. 

6.2 Correlation length and osmotic pressure of semi-dilute polyelectrolyte solutions 

In a dilute solution, polyelectrolyte chains have a rodlike conformation at a local scale 
smaller than the persistence length £ p and the persistence length is always larger (or of the 
same order of magnitude) than the screening length of the electrostatic interactions. In a 
semidilute solution, the polymer chains overlap and, if we assume that they do not form 
ordered phases, they form a temporary network with a mesh size £. In the absence of salt, 
the chains also have a rodlike conformation at the scale £ and their persistence length is, 
as shown below, larger than the mesh size. The mesh size can then be calculated from 
a scaling argument by imposing that it is equal to the isolated chain radius (2.4) at the 
overlap concentration c*. (For simplicity, we suppose here that locally the polymer shows 
local gaussian statistics, the other cases can be treated in a similar way) . Equivalently, a 
pure geometrical argument can be used, by imposing a close packing condition for chain 
subunits of size £, c = (C/Ce)d/^- I* 1 an y case, £ is given by [14,39] 




(6.1) 
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The number of monomers within a volume of size £ , or correlation blob, is then = 

For weakly charged polyelectrolytes, the counter ions are free and behave roughly as an 
ideal gas. The screening length due to the counterions is then k^ 1 = (47r£ s /c) -1 / 2 . This 
screening length is of the order of the mesh size for strongly charged polyelectrolytes (/ ~ 1) 
but it is larger than £ for weakly charged polyelectrolytes. The polymer itself, however, 
contributes to the screening of the electrostatic interactions and the actual screening length 
is smaller than the counterion screening length. From a purely electrostatic point of view 
each cell of size £ containing monomers is neutral and can be considered as independent 
of all the others. It is thus reasonable to assume that the effective screening length is of 
the order of the mesh size £ [43,75]. 4 

As salt is added to the solution, the screening length decreases and becomes dominated 
by salt (k 2 = 8nn£ B ) when the Debye-Hiickel screening length of the salt k^ 1 is smaller 
than the mesh size £. When the screening is dominated by the salt, the electrostatic 
persistence length of the chains decreases with the salt concentration. As long as the 
persistence length is larger than the mesh size £, the mesh size is still given by equation 
(6.1). If the persistence length is smaller than the mesh size, the polyelectrolyte solution 
behaves as a neutral semi-flexible polymer solution and the mesh size or correlation length 
can be calculated from the radius in a dilute solution using scaling arguments. If the 
Odijk-Skolnick-Fixman statistics is assumed [76,43] 

£ ~ 6£ eC - 3 /V/ 4 &- 3 / 2 , (6.2) 
while if the persistence length is assumed equal to the screening length [14,39] 

£ ~ £ e 3 / 4 c- 3 /V/2 6 -3/2 (6.3) 

This is somewhat similar to the assumption that, in neutral polymers, the hydrodynamic 
screening length is equal to the static correlation length. 

42 



Experimentally, the important feature of semidilute polyelectrolyte solutions is the fact 
that, if the salt concentration is low enough, the structure factor S(q) has a peak at a finite 
wavevector q*. At small wavevectors, the structure factor, because of the electroneutrality 
of the solution, is dominated by the small ions, that behave as an ideal gas. In the 
absence of salt, S(q = 0) = 1/ f. If the salt concentration n is larger than the counterion 
concentration, S(q = 0) = 2n/(/ 2 c). At large wavevectors, corresponding to distances 
smaller than the mesh size £, the polyelectrolyte has a rodlike behavior and the structure 
factor is given by S(q) = £, e /(qb 2 ). When the wavelength g _1 is equal to the mesh size, 
the value of the structure factor is S* = = c . It can be checked explicitly 

that this value is larger than the thermodynamic value S(q = 0) and the structure factor 
must thus have a peak at a finite wavevector q*. The existence of the peak is therefore 
related to the very small compressibility of the small ion gas which gives S(q = 0) due to 
the electroneutrality constraint. This discussion is valid as long as the counterions and the 
salt ions are uniformly distributed throughout the solution, i.e. as long as the screening 
length is larger than the mesh size of the solution. In the salt dominated regime, where 
the salt screening length is smaller than the mesh size of the solution, the counterions are 
confined in a sheath of size around each polymer and the solution essentially behaves 
as a neutral polymer solution where the structure factor decays monotonically. The peak 
thus disappears in the salt dominated regime. At the crossover between the two regimes, 
the q = value of the structure factor increases very sharply. When the effect of salt is 
not dominant, the peak position q* defines the correlation length of the solution and one 
expects that it is of the order of the inverse of the mesh size q* ~ l/£. The position of 
the peak of the structure factor increases thus as c 1 / 2 in a semidilute solution (and as c 1 / 3 
in a dilute solution as explained above). It is also important to note that the correlation 
length £ and thus the peak position are roughly independent of the salt concentration as 
long as the Debye-Hiickel screening length is smaller than £. These results are in rather 
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good agreement with experiments. 

In the osmotic regime where the Debye-Hiickel screening length is smaller than the mesh 
size, the small ions are uniformly distributed and the osmotic pressure of the solution is 
the same as that of a gel as discussed in section 5, it is dominated by the counterions. It is 
equal to tv = kgTfc if the counterion concentration is smaller than the salt concentration, 
(fc > n) and to 7r = /c^T/ 2 c 2 /(4n) if fc < n. In the salt dominated regime k~ 1 < £, the 
polyelectrolyte chains behave as neutral semiflexible chains with a persistence length £ p 
and an effective diameter equal to the screening length (up to logarithmic corrections). 
The excluded volume between Kuhn segments of size £ p is of order £ 2 k _1 . If the persistence 
length is larger than the mesh size (or of the same order), the osmotic pressure varies as 

7v/k B T = c 2 K-^- 2 b 4 (6.4) 

At the crossover between the osmotic and the salt dominated regime, the osmotic pressure 
varies very rapidly and the crossover is not smooth. When the persistence length becomes 
small enough compared to the mesh size, the excluded volume correlations are relevant and 
the osmotic pressure increases as n ~ c 9 / 4 . At very high ionic strength, the electrostatic 
interaction is very weak and cannot be approximated by an excluded volume between Kuhn 
segments. In this weak coupling regime there no longer is an electrostatic persistence length 
(it is shorter than the screening length) and the polyelectrolyte chains should be considered 
as flexible chains where the monomers interact via the short range Debye-Hiickel potential 
with an effective excluded volume v e \ = / 2 /2n. The smooth crossover between the salt 
dominated and the weak coupling regimes occur when K^ e ~ 1 i.e. when the screening 
length becomes equal to the electrostatic blob size. 

6.3 Electrostatic rigidity in semidilute solutions 

In this section we discuss the conformation of polyelectrolyte chains in a semidilute 
solution and the effect of the interactions between different chains on the persistence length 
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of a polyelectrolyte. For simplicity we consider only a strongly charged polyelectrolyte for 
which the bare persistence length £ and the distance between charges A are such that 
A 2 > £ £ B so that the Odijk-Skolnick-Fixman theory can be applied in a dilute solution. 
We also assume that the bare persistence length is much smaller than the mesh size of the 
semidilute solution. The results can easily be extended to weakly charged polyelectrolytes 
if necessary by renormalising the bare persistence length to the blob size. 

Two types of theories have been proposed to describe chain conformation in a semidi- 
lute solution. The early work of Odijk [76,77] makes the assumption that the interchain 
interactions have a negligible influence on the persistence length so that the persistence 
length is given by the usual Odijk-Skolnick-Fixman theory with the relevant screening 
length. A scaling theory for the chain conformation has been built on this assumption. It 
was however later pointed out by Witten and Pincus [78] that if a polyelectrolyte solution 
is isotropic, the electrostatic interaction energy between two chains can be reduced by a 
bending of the chains to avoid each other and thus that the interaction energy between 
different chains reduces the persistence length. A persistence length equal to the mesh size 
of the solution is obtained in the osmotic regime where the Debye-Hiickel screening length 
of the salt is larger than the solution mesh size. In a revised version of this theory, in the 
salt dominated regime, the corrections to the persistence length due to the interactions 
between chains are found to be small. 

The simplest approach [79] to the persistence length of a polymer chain in a semidilute 
solution is to use linear response theory. Within the framework of linear response, an 
effective pair interaction between monomers can be introduced. The Fourier transform of 
the effective interaction is related to the structure factor S(q) of the solution by v e fj(q) = 
vdh(q) (1 — c S(q)vDH(q)). Following the lines of Odijk argument, the effective persistence 
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length can then be calculated as 




ftpg (g = °) 
16ivk B TA 2 



(l-cS(q = 0)v DH (q = 0)). 



(6.5) 



This relation gives the persistence length as a function of the q = value of the struc- 
ture factor which is a thermodynamic quantity. If the explicit form of the Debye-Hiickel 
potential is substituted, the persistence length is l v = £ + £ OSF (l — 4n£ B cS(q = 0)/k 2 ). 
A reduction of the persistence length from the single chain value induced by interchain 
interactions is therefore predicted. However this result can be used only in the limit of 
linear response, i.e. if the interaction between chains is smaller than kgT. The interaction 
between two rods crossing at perpendicular angle is 



The linear response result is thus valid if (3 « 1 which corresponds to the weak coupling 
regime for rigid polyelectrolytes. 

No direct calculation of the effect of interchains interactions on the persistence length 
seems available in the strong coupling regime j3 » 1. Some insight into the problem 
can however be gained by looking at the following simpler two dimensional problem: one 
semiflexible test chain lies in a plane, and interacts with all the other chains replaced by 
infinite rods perpendicular to this plane with a concentration per unit area V. It is a 
good approximation to replace the other chains by infinite rods as long as their persistence 
length is larger than the screening length. In the limit where the Debye screening length 
is smaller than the distance between the chains, the solution can be considered as a quasi 
ideal gas of rods. When the test chain bends, the distribution of the rods around the test 
chain is changed and this contributes to the bending energy. The persistence length can 
be calculated following the lines of the work of Odijk in a dilute solution, the result is 



OSF 



(6.7) 



46 



where the function f((3) is defined as 



f(P) 




dzexp(-/3exp(-z)) -(3z 2 exp(-2z) + exp(-z)(-z 2 + z + 1) (6.8) 



The correction to the Odijk value is always negative which corresponds also to a reduction 
of the persistence length from the isolated chain value. In the weak coupling limit, f{(3) = 
3/9 2 /8 in agreement with the linear response theory. In the strong coupling limit, the 
interaction is very similar to a hard sphere interaction and f((3) = (ln/3) 2 varies only 
logarithmically with the interaction strength. The correction to the single chain persistence 
length due to interchain interactions is small. 

These results obtained in two dimensions find a simple interpretation in terms of indi- 
vidual deflections of the test chain by the perpendicular rod. When the test chain interacts 
with one rod at a distance r, it bends and its orientation changes over a distance by 
a finite angle that can be shown to be equal to 



The persistence length can then be obtained by summing the individual deflections weighed 
by the Boltzmann factor related to the direct electrostatic interactions between the chains. 
The persistence length given by equation (6.7) is then found up to a numerical prefactor. 
The independent deflection model is easily generalisable to real three dimensional chains. 
The chains do not always cross at a finite angle, but the conformations where two chains 
are almost parallel are energetically costly, since the electrostatic interaction is very high 
in this case (it increases linearly with chain length). This leads to a persistence length 



where a is a numerical constant of order unity. The important result is that the correction 
is small and thus that in the salt dominated regime, the reduction of the persistence length 
due to the interactions between different chains can be neglected. 



6 ~ /trexp(— K,r). 



(6.9) 




(6.10) 
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The situation is very different in the osmotic regime, the chains interact very strongly 
and the distance between interacting chains is always of the order of the mesh size £ ~ 
the screening length k 1 is also of order £. There is no precise theory of the 
persistence length in the osmotic regime, the two following qualitative arguments suggest 
however that the persistence length is of the order of the mesh size £. As the distance 
r between interacting chains is of order the deflection angle given by equation (6.9) 
is always of order one. Whenever two chains cross they bend by an angle of order unity 
and loose the memory of their orientation. The persistence length is thus of the order of 
the distance between crossings i.e. of the order of the mesh size £. The second argument 
is based on the strong structuration of the solution in the osmotic regime. The structure 
factor has a strong peak and with a very rough approximation the chains can be thought of 
as lying on a lattice with a lattice constant of order £. The total electrostatic energy of this 
lattice does not depend of the bending of the chains over distances larger than £. When 
the chains are bent on a scale £, entropy is gained but at the expense of bending energy 
associated with the bare persistence length £ . The entropy is clearly dominant when £ < £ 
which leads to a persistence length of order £. An important consequence is that in saltfree 
semidilute solutions, both flexible and rigid polyelectrolytes have a persistence length of 
the order of the correlation length. It is also important to notice that the persistence 
length of rigid polyelectrolytes is predicted to vary non-monotonically with ionic strength. 
In the osmotic regime, it is dominated by interchain interactions and increases with ionic 
strength. In the salt-dominated regime, it is dominated by intrachain interactions and 
decreases with ionic strength. There is no clearcut experimental evidence of this non 
monotonic variation of the persistence length with ionic strength [80] . 

6.4 Concentrated solutions of flexible polyelectrolytes 

The semidilute regime discussed in the two previous subsections corresponds to chains 
which are rodlike at intermediate length scales and that can be described locally by the 
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electrostatic blob model of section 2. At a higher concentration however, the mesh size 
of the solution is smaller than the electrostatic blob size. Electrostatic interactions only 
play a minor role in this regime and the polyelectrolyte chains remain gaussian at all 
length scales. The electrostatic blob size £ e is equal to the correlation length £ when the 
electrostatic blobs are close packed. If the local structure of the blobs is gaussian, this 
occurs at a concentration 

c « = 6 -3/2/3 (<fl/6) l/3 ( 6 .11) 

In the concentrated regime c > c**, the elecrostatic interactions are weak and the correla- 
tions in the solution can be studied within the framework of the so-called random phase 
approximation, as first suggested by Borue and Erukhimovich [81,82]. This leads to a 
structure factor (monomer-monomer concentration correlation function) 

1 1 2 4tt£ b / 2 , C10 . 

S(q) S (q) q 2 + K 2 

Here So(q) is the structure factor of gaussian chains, given by the Debye function, and in 

the following it is written as So(q)~ 1 = (iVc) _1 [l + g(Nq 2 b 2 /6)]. The second and third 

virial coefficients that describe the interactions between neutral monomers (i.e. that would 

describe the interactions for / = 0) are, respectively, v and w. The second virial coefficient 

v is the excluded volume v and vanishes at the 9 temperature. The third virial coefficient 

w 2 is assumed here to be positive. The last term comes from the screened electrostatic 

interactions. The screening is due to all the small ions and the Debye-Hiickel screening 

length is defined here by k 2 = 4n£ B (fc + 2n). The statistics of the electrostatic blobs 

remain gaussian if the excluded volume is small enough, namely if the electrostatic blob 

size £ e is smaller than the so-called thermal blob size & ~ b^/v. 

Within the RPA approximation (6.12), the inverse osmotic compressibility of the neu- 
tral solution is l/(Nc) + v + w 2 c. That of the charged solution is l/(Nc) + v + w 2 c + 
47r£ B f 2 /K 2 . Hence, in the small wavevector limit, the electrostatic interactions give rise to 
an additional excluded volume v e \ = Anl B f 2 / k 2 . 
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At high wavevectors, the chains are gaussian and the structure factor decays as S(q) ~ 
12/ (q 2 b 2 ). At small salt concentrations, the structure factor given by (6.12) has a peak at 
a wave vector q* given by 

9* 2 + « 22 = — ( 6 - 13 ) 

where </ is the derivative of g with respect to (Nq 2 b 2 /6), equal to 1/2 when (Nq* 2 b 2 /6) is 
large and to 1/3 when it is small. The peak position is independent of the solvent quality. 
In the absence of salt, it scales as q* ~ / 1 / 2 c 1 / 4 and crosses over smoothly to the peak 
position in the semidilute range q* ~ l/£ given by equation (6.1) at the concentration 
c**. As the salt concentration is increased, the peak shifts towards zero wavevector and 
disappears when > (24n£ B f 2 c)/(g f b 2 ). At higher ionic strength the polymer solution 
behaves as a neutral polymer solution and the structure factor decreases monotonically 
with the wavevector. The structure factor given by equation (6.12) is in good agreement 
with neutron scattering experiments on weakly charged polyacrylic acid or polymethacrylic 
acid solutions in water 



The osmotic pressure of the solution can also be calculated from the RPA (one loop) 
approximation. In the absence of salt it is equal to 

— — — = fc -\- -vc 2 + -W 2 C 3 - 4f 3 /2 c 3/4 £ 3/4 6 -3/2 (6 . 14) 
tzj^T 2 3 

where A is a numerical constant of order unity [81]. The first term is the osmotic pressure 
of the counterions discussed extensively above; the two following terms are due to the 
non electrostatic interactions between monomers. The last term is due to the charge 
fluctuations and is very similar to the so-called polarisation pressure of simple electrolytes. 

Finally, the cooperative diffusion constant of the solution can be calculated from the 
structure factor using linear response theory and describing the hydrodynamic interactions 
by the so-called Oseen tensor [84] . When the electrostatic interactions dominate over the 
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excluded volume interactions, the cooperative diffusion constant is given by 

D = + „ 2 /« 2 )(2 + K */a 2 )- 1/2 (6-15) 

OUT] K l 

where r\ is the solvent viscosity and where the wave vector a is defined as a = (g* 2 + ^ 2 ) 2 . 
In a first approximation, the cooperative diffusion constant varies as D ~ / 3 / 2 c 3 / 4 (/c + 
2n)~ 1 . It is thus a non monotonic function of the monomer concentration c. It increases at 
low concentration, reaches a maximum when the counterion concentration is of the order 
of the salt concentration and decreases at higher concentration. The unusual decrease 
at high concentration is due to the coupling to the counterion gas which dominates the 
compressibility of the solution. This might explain the unusual variation of the diffusion 
constant observed on polyacrylic acid solutions and gels. 

A summary of the various regimes expected for polyelectrolyte solutions in the semidi- 
lute and concentrated ranges is given in figure 2. 

6.5 Mesophase formation in poor solvents 

In a poor solvent below the temperature, the polyelectrolyte chains are subject to 
two antagonist forces, the attractive (Van der Waals like) interaction and the repulsive 
Coulombic interaction. As explained for isolated chains in section 2, the attractive inter- 
action is dominant at short length scales and induces a local collapse of the solution. If 
the salt concentration is not too high, phase separation is not expected at a macroscopic 
scale but at a mesoscopic scale. The polymer solution is expected to form mesophases 
where the polyelectrolyte concentration is not homogeneous and where polymer dense and 
polymer dilute regions are periodically arranged. These phases are very similar to the 
mesophases observed in surfactant or block copolymer systems and various symmetries are 
possible (lamellar, cubic and hexagonal phases) [85,86]. Qualitatively, the mesophases are 
stabilised by the translational entropy of the counterions. If a macroscopic phase separa- 
tion occurs, the separated phases must be electrically neutral and the translational entropy 
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of the counterions is low. When a mesophase is formed, the electroneutrality can be locally 
violated, this is associated to a cost in electrostatic energy but the translational entropy 
of the counterions is tremendously increased compared to a neutral dense state. 

In the vicinity of the 9 temperature, mesophase formation can be investigated using 
the same random phase approximation as in the previous section. The structure factor is 
given by equation (6.12) where the excluded volume v = — 6 3 r is now negative, r being 
the temperature shift from the 9 point. The structure factor has a strong peak at a 
wavevector q* and this peak diverges as the temperature is lowered (v becoming more 
negative) indicating the formation of a periodic structure with a period 2% / 'q* . In the 
absence of salt, the spinodal for mesophase formation obtained from the divergence of the 
peak of the structure factor occurs when 

-+vc + w c +-g(Nq b /6) + — gf p J =0 &.lb) 

The qualitative shape of this spinodal line in a temperature r-concentration c plane is 
shown in figure 3. The temperature and the concentration at the maximum of the spinodal 
line are given by r m ~ / 2 ^ 3 (^ B /^) 1 ^ 3 and c m ~ 

6 -3^2/3^ b / 6 )1/3_ Upon addition of salt, 
the period of the mesophase increases and when q* = (k a = (24n£ B f 2 c)/(g'b 2 )), the 
mesophase transition becomes a macroscopic phase separation very similar to the demixing 
transition of a polymer in a poor solvent. 

In this weak segregation limit the symmetry of the phases and the topology of the 
phase diagram have been calculated by the Moscow group [87] by studying the non linear 
fluctuations of the solution. 

At a temperature much smaller than the 9 temperature, in the strong segregation 
regime, the interface between the dilute and dense polymer regions is sharp and can be 
described in terms of a surface tension. In a very dilute solution micelles are expected to 
form and their characteristics have been studied in references [82,88]. At higher concentra- 
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tions periodic phases are predicted. A partial phase diagram considering only the lamellar 
ordering has been determined by Nyrkova et al. [89] by a direct minimisation of the mean 
field free energy. 

Most experimental results seem to be obtained outside the region of the phase diagram 
where mesophases form. To our knowledge, mesophase formation has been reported only 
in one case [90] and no detailed study of the mesophase structure seems to exist. 

7. Dynamical properties 

In this section, we briefly discuss some dynamic properties of polyelectrolyte solutions. 
Only two quantities are considered: the mobility and the electrophoretic mobility of a chain 
in a very dilute solution and the viscosity of the solution. In an electrophoresis experiment, 
the hydrodynamic interactions are screened and the electrophoretic mobility is independent 
of the molecular weight. The separation of molecules of different size (such as DNA) is thus 
not possible by simple electrophoresis. The concentration variation of the viscosity of a 
polyelectrolyte solution is very different from that of an organic neutral polymer solution. 
This effect has been known experimentally for a long time (it is sometimes referred to as 
the polyelectrolyte effect) but it remains poorly understood. 

7.1 Mobility and electrophoretic mobility of a single charged chain 

The theory of transport properties in dilute electrolyte solution has been worked 
out long ago by Debye, Onsager and Falkenhagen (see [91], [92], and historical references 
therein). The existence of long range Coulomb interactions between the species in solution 
results in a deviation of the transport coefficients from the pure solvent values (ionic mobil- 
ity, conductivity, electrophoretic mobility and viscosity) that increases as the square root 
of the concentration, rather than the linear increase which is the rule for neutral species. 
The origin of this anomalous behaviour can be traced back to two effects known as the 
relaxation effect and the electrophoretic effects. Although these effects have different phys- 
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ical origins, they both produce a variation of the transport coefficients as the square root 
of the ionic concentration. 

The relaxation effect can be described as follows. A test charge placed in an ionic 
solution is surrounded by a spherically symmetric polarisation cloud where charges of 
oppposite sign are predominant, that can be described by the standard Debye-Hiickel 
theory. If the test charge moves at a constant velocity U with respect to the solvent, the 
polarisation cloud is slightly distorted, and the charge experiences an additional electric 
field that tends to slow down its motion. The resulting force therefore appears as an 
electrostatic friction on the test charge. More generally, the moving charge Q (position 
R, velocity U) in an electrolyte solution can be shown (see Appendix B) to produce an 
electric field inside the solution of the form 



The ionic mobility (q is f° r simplicity assumed to be identical for all small ions. 

The electrophoretic effect is important only for the description of ionic motion driven 
by an external electric field (electrophoresis). In that case, the force that drives a charge 
motion also acts on the surrounding polarisation cloud. The solvent velocity field that 
results from this body force acting on the polarisation cloud modifies the relative velocity 
of the charge with respect to the solvent creating an increased dissipation. In other words, 
the volume of solvent within a distance k~ 1 of the moving charge is dragged along as the 
charge moves, with a friction coefficient of order r]K~^ , where rj is the solvent viscosity. In 
the language used to study polymer dynamics, the electrophoretic effect can be described 
as a screening of the hydrodynamic interactions when the motion is created by an electric 



E(r) 



T(r-R)U 



(7.1) 



where the a (3 component of the matrix T is 




(7.2) 
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field. The hydrodynamic interactions are usually described in terms of the Oseen tensor 
H(r) [6], which gives the velocity field that results from the application of a local force 
Fo5(r) localized at the origin. When the force results from an electric field Eq applied to a 
charge Q located at the origin, it is easily shown (see Appendix B) that, if the polarisation 
cloud is described within the Debye-Hiickel approximation, the velocity field is given by 

exp(-Kr)H(r)QE (7.3) 

This velocity field is screened in a similar way as the electrostatic interaction and decays 
as exp(— Kr)/r. The hydrodynamic interaction with neighboring particles is also screened 
over a distance 

The extension of these concepts to macroions has been the subject of numerous con- 
tributions. Much of this work however is concerned with spherical ([24], and references 
therein) or rodlike [93] particles. The case of flexible polymers has been considered by 
Hermans [94] Manning [92] and Muthukumar [95]. In reference [92], the relaxation and 
electrophoretic effects are treated on a different footing. In reference [95], the relaxation 
effect is ignored. In the following, we present a simplified description that allows for a 
simultaneous treatment of the two effects at a microscopic level. Our conclusions are es- 
sentially similar ^ of [92]. They differ from those of reference [95], which we believe to be 
in error. 

We consider here two types of experiments. In a sedimentation experiment, the external 

force is the same on all monomers and has no effects on the small ions. The chain mobility 

H relates the velocity of the center of mass of the molecule vq to the total force acting 

on the molecule F: vq = fiF. In the electrophoresis experiment, the external force is due 

The result that is arrived at in [92] can be obtained within our analysis by ignoring the sec- 
ond term in the r.h.s. of equation (7.6). The result that is obtained for the electrophoretic 
mobility is then (i e = fe(H s o) /(l + fe(H s o)To), which is only slightly different from (7.9), 
since (H s q) is independent of molecular weight. A problem would arise, however, if the 
same approximation were made in the calculation of the mobility. 
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to the applied electric field E acting both on the charged monomers and on the free small 
ions. The electrophoretic mobility is defined by vq = /U e E. 

The model that we use to describe the dynamics of the charged polymer is a simple 
extension of the standard Zimm model of neutral polymers [6]. The velocity and position 
of monomer i are denoted respectively by R^ and Vj, the velocity field of the solvent is 
u(r). The equation of motion of monomer i, bearing a charge fe results from a balance 
between a friction force and all the other forces acting on the monomer: 

C(vi-u(R z ))=F z + /eE r (R 2 ) (7.4) 

The monomer solvent friction coefficient is (. The force Fj is the sum of an external force 
F ex t (the sedimentation force or the applied electric field), a Langevin random force 6>j and 
of an intramolecular force F po i^. The last term E r is the electric relaxation field. 

The velocity field of the solvent is given by an Oseen-type relation 

u (R-i) = ^2 _ Rj)i s (Ri - Rj)F ex t + Fpoij + 9j] + a s (Ri - -Rj)/eE r (Rj)] ^ ^ 

3 

the screening factor s(r) is equal to exp(— nr) in the case of electrophoresis and to unity 
otherwise. If we assume that the solvent velocity field varies slowly over a scale k -1 , the 
relaxation field E r (Tti) on monomer i is obtained from (7.1) as 

E r (R i ) = -J]T(R i -R i )(v j -u(R j )) • (76) 

3 

These equations can be treated in the spirit of the Kirkwood-Rieseman preaveraging ap- 
proximation, by replacing the tensors H(Rj — Rj), s(Rj — Rj)H(R — R ? -) and T(R — Hj) 
by their equilibrium average values. Introducing 

H =N- 1 <]TH(R Z -R 7 )> 

ij 

H s0 =iV- 1 < ^exp(- K |R z - R J |)H(R Z - Rj) > 



T =iV- 1 <^T( J R z - J R J )> 

ij 
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(7.7) 



the mobility /i and electrophoretic mobility fi e are expressed as: 



V = ^ T \Ho-H s0 fe^(l + fe^)- 1 



(7.8) 



lie = feH s0 (l + fe— ) 



(7.9) 



The mobility is given by the usual Kirkwood-Riseman result [6], with a correction that 
arises from the relaxation effect. This relaxation correction reduces the mobility. Using 
equation (7.2), Tq can be computed as 



where S(q) is the structure factor of the chain (normalized so that S(q = 0) = N). From the 
second line of (7.10) it is seen that the effect of the relaxation field is independent of molecu- 
lar weight whatever the chain conformation and is small when the charge fraction / is small. 
For example, if the polymer is rodlike over a distance /eTg/C = / 2 ^ B / (3&)(Co/0- The 
averaged screened Oseen tensor, H s q, is also independent of molecular weight and propor- 
tional to — ln(«6). The mobility is thus dominated by the Kirkwood Riseman contribution, 
and is essentially the same as for a neutral chain with the same conformation. The re- 
laxation correction becomes important for strongly charged rodlike chains in the vicinity 
of Manning condensation threshold. In this case, the mobility is reduced by a factor of 
order unity from the Kirkwood Riseman value. The approximations that we have made 
are only valid for long chains, R» In view of the complex form of equation (7.8), it 
is not surprising that experimental results for polyelectrolyte diffusion constants are often 
intermediate between the "non-draining" and "free draining" limits [96]. 

The electrophoretic mobility fi e , on the other hand, is in any case independent of molec- 
ular weight because of the screening of hydrodynamic interactions. For weakly charged 




(7.10) 



frf^N-l < ]Texp(- K |R, - Rj\) > 
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polyelectrolytes, the relaxation correction is small and the electrophoretic mobility can be 
rewritten in terms of the structure factor of the chains as 

^3&/ d * SW (FT^) < 7 - n > 

For a rodlike molecule, the electrophoretic mobility decreases logarithmically as the ionic 
strength is increased fj, e = — 3^25 ln/tfe. For a flexible chain, using the stretched chain 
model of section 2 to describe the local conformation of the chains, we obtain \i e = 
— (fb/£ B ) 1 ^(e/r]b)\iiK$, e , the increase with the charge fraction is weaker and the elec- 
trophoretic mobility also decreases logarithmically with ionic strength. The general fea- 
tures of these theoretical predictions are in reasonable agreement with experimental ob- 
servations. A detailed comparison between theoretical and experimental results is given in 
[92]. 

7.2 Viscosity of poly electrolyte solutions 

The viscosity of polyelectrolyte solutions has been known very early to exhibit a be- 
haviour qualitatively different from that of neutral polymer solutions (see [1] , and references 
therein). This behaviour was first described by Fuoss [97] using an empirical law, which 
gives the reduced viscosity r\ r = (77(c) —r) s )/rj s c as a function of the polymer concentration 

c 

* = TTlU ' (7 ' 12) 

Detailed studies [98] show that the constant A, obtained by extrapolating the results to 
vanishing concentration, is proportional to iV 2 / 2 , while the ratio A/B is proportional to 
Nf. 

It was later realized that this law is satisfactory only for high enough concentrations 
and in saltfree solutions, in which case it reduces to 7] r ~ c~ 1//2 . At low concentration, 
the behaviour of r\ r is actually nonmonotonic [99]. The reduced viscosity r] r presents a 
peak at a finite concentration, of the order of the salt concentration. This peak disappears 
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when enough salt is added to the solution. A phenomenological law, recently proposed 
by Cohen, Priel and Rabin [100], seems to account well for the experimental observations. 
According to these authors, the intrinsic viscosity is given by 

Vr = ^ (7-13) 

where k is the Debye screening parameter, and A\ is proportional to the molecular weight. 

Finally, it is important to note that this anomalous behaviour of polyelectrolyte solu- 
tions can be rationalized at low concentrations by using the method of isoionic dilution 
[1,101]. In this method, the viscosity is studied as a function of polymer concentration for 
a given value of the ionic strength (i.e. as the solution is diluted, salt is added in order 
to keep the ionic strength constant). The increase of the reduced viscosity at low c is 
then linear, as for neutral polymers. In neutral polymer solutions, this linear behaviour is 
usually represented in the form r\ r = [77] + &#[??] 2 c, where [rj] is the intrinsic viscosity and 
kfj is the Huggins constant. In polyelectrolyte solutions, [rj\ is large and varies strongly 
with the ionic strength /. The data of [101], for exemple, can be represented in the form 
W = [v](I = °°) + CT -1 / 2 [102]. The slope of the linear term, kjj[rf\ 2 , is also much larger 
than the corresponding quantity in neutral polymers. The Huggins constant kjj varies as 
J -1 / 2 [102] 6 . This indicates the importance of interactions between chains, and also shows 
that, as c is increased a crossover to the "interaction dominated" behaviour described by 
(7.13) rapidly takes place. 

In short, three different concentration regimes can be distinguished in the behaviour 
of the reduced viscosity as a function of c (figure 4) 7 . In the low concentration region 
(left side of the peak), i] r is an increasing function of c. This region is most usefully 

In fact, the Huggins constant is not a particularly useful concept in a system where inter- 
actions that are not of hydrodynamic origin are expected to be important. 

This distinction is a simplification of the description proposed by Wolff [102], and seems 
valid at all but the highest concentrations. 
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characterized by isoionic dilution curves. The peak region corresponds to comparable salt 
and counterions concentrations (fc ~ n). Finally, in the more concentrated region, rj r is a 
decreasing function of c, r\ r ~ c~ v with y < 1. Obviously, this description is valid only in 
solutions where the salt concentration is small enough, since if n is large the behaviour of 
usual neutral polymers must be recovered. 

At present, a comprehensive theory of polyelectrolyte viscosity is lacking. A number of 
interpretations, however, have been put forward to explain the unusual behaviour of these 
solutions. We will distinguish here two types of models based either on the conformation 
of the polymers or on the interactions between different chains. 

In conformation based models [14,103,104,105] (which have been formulated for flexi- 
ble polyelectrolytes), it is assumed that the viscosity of the polyelectrolyte solution is only 
a consequence of the structural properties of the solution. The calculation proceeds by 
analogy with the well understood case of neutral polymers. For semidilute solutions, the 
viscosity can be computed by applying the Rouse or reptation theories to a chain of corre- 
lation blobs of size £, containing G^ monomers. The correlation length of polyelectrolyte 
solutions has a variation with concentration very different from that of neutral polymer 
solutions, so that the behaviour of the viscosity is also qualitatively different. This idea was 
first applied by de Gennes and coworkers [14]. These authors interpreted the [77] ~ A^c~ 1//2 
behaviour in saltfree semidilute solutions given by Fuoss law as corresponding to Rouse 
dynamics for a chain of correlation blobs of size £ ~ k (equation (6.1)) with G^ = g(£/£e) 
monomers. The viscosity is then given by 



This analysis was extended by Dobrynin et al to the case where salt is added to the solution 




(7.14) 
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[75]. In that case, £ is given by (6.3), = c£ 3 and the viscosity is 



= A^c 5 / V 3 / V 7 V/ 2 . 



(7.15) 



The crossover between (7.14) and (7.15) takes place when the mesh size is equal to k , 
like the crossover between (6.1) and (6.3). Note that equation (7.15) predicts a maximum 
in the reduced viscosity rj r (c) for fc = n. 

Rubinstein and coworkers [75] also give an empirical justification for the observation 
that semidilute polyelectrolyte solutions obey Rouse dynamics over a broader range of 
concentrations above the overlap concentration than neutral polymer solutions. In neutral 
polymers, it is observed that entanglements become important when each chain inter- 
acts with about 10 other chains [106]. If this criterion is applied to polyelectrolytes, the 
concentration c rep t for the onset of reptation behaviour is given by 



This can be rewritten as c rep t^(c rep t) — N/100. The overlap concentration c*, on the 
other hand, is given by c*£(c*) 3 ~ N. In polyelectrolyte solutions, c£ 3 (c) decreases much 
more slowly with increasing concentration c than in neutral polymer solutions, so that the 
onset of reptation is expected at higher concentrations. If equation (6.1) is used, the ratio 



Another useful way of computing the viscosity of neutral polymer solutions, which 
allows a description of the dilute/semidilute crossover, is to write it in a scaling form 



where rj s is the solvent viscosity, F is a scaling function and c* is the overlap concentration. 
The function F can be expanded as a power series (F(x) = 1 + a\x + ci2X 2 ..) for small 
values of its argument, i.e. in dilute solutions. In the semi-dilute regime, the functional 




(7.16) 



Crept /c* is of order 10 , while for neutral polymers it is of order 10 . 




(7.17) 
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form of F depends on the mechanism (Rouse or reptation) that governs stress relaxation. 
In this formulation, the specificity of polyelectrolyte solutions is related to the fact that 
the interactions between monomers depend on the polymer concentration. Therefore c* in 
(7.17) must now be interpreted as the overlap concentration corresponding to a polymer 
having the same interactions as the charged polymer at concentration c. For exemple, 
when the mesh size is larger than we can write 

c*(c) = N/R(cf. (7.18) 

Where R(c) = K _1 (f £e«) 3/5 , 

given by equation (4.12), is the radius of an isolated chain 
with this interaction. The viscosity is obtained by writing rj = rj s F(c/c*(c)), and assuming 
Rouse behaviour r\ ~ N. This gives 

/ n \ 5 / 4 

V = Vs (-^y) = (7.19) 

which is identical to (7.15). 

Another " conformational" interpretation of polyelectrolyte viscosity was given by Reed 
[104], using similar arguments. His calculation assumes that c is in the vicinity of c*(c), 
so that the x 1 contribution in the function F(x) dominates. The formula used for c* is 
appropriate for short chains or large It corresponds to an ideal behaviour for a chain 
with a persistence length k~ 1 

R(cf = 2(N/g)i e n- 1 . (7.20) 

This leads to a viscosity 

r, ~ r, s (c/c*f ~ Vs N^ e /gf. ( 7 - 21 ) 

(Note that the same result would be obtained by requiring Rouse behaviour for the viscos- 
ity, i.e. F(c/c*) ~ N) Again, the behaviour r\ r ~ /A/"c -1 / 2 is obtained in the absence of 
salt. With added salt, the experimentally observed behaviour (7.13) is reproduced. The 
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main difference between the interpretations of Reed and of Dobrynin et al. is the expres- 
sion used for R(c). The crossover from (7.21) to (7.15) can be expected to take place as 
c increases and R(c) crosses over from the ideal expression (7.20) to the Flory expression 
(2.5). 

Finally, we note that in all conformation based models the c -1 / 2 behaviour of r] r , 
corresponds to situations where the role of added salt is negligible. It is interesting that 
the value of the constant A in Fuoss law (7.12), which corresponds to a (generally incorrect) 
extrapolation of this behaviour to c = 0, yields a value A ~ N/R? ~ iV 2 / 2 that corresponds 
to the intrinsic viscosity for a saltfree, dilute solution of polyelectrolytes (R ~ iV/ 2 / 3 ). This 
coincidence, however, is not explained theoretically. At high concentrations, entanglements 
become important, and the reduced viscosity increases with concentration [102,75] 

The conformation based models appear to provide a qualitatively correct description 
of flexible polyelectrolyte viscosity (see [104,75] for a comparison with experimental data). 
They only apply, however, to sufficiently concentrated solutions, near or above overlap 
concentration. The "polyelectrolyte behaviour" described by equations (7.12) or (7.13) 
is also observed in solutions of short chains, or even for charged spherical particles [107]. 
In [108], it was suggested that this behaviour could be described on the basis of the liq- 
uid state "mode-coupling" theories developped for the study of spherical charged colloids 
[109]. Unfortunately, these theories are analytically tractable only in the limit of weak 
electrostatic coupling between the particles, which is certainly not realized in dilute poly- 
electrolyte solutions (see section 6.1). In this weak coupling regime, it was shown that 
the viscosity of a suspension of spherical charged particles (charge Z, concentration cq, 
hydrodynamic radius R) is given by 

rj ~ rj s Rc 2 Z A £ 2 B K- 3 (7.22) 

The authors of reference [108] noticed that this behaviour is similar to the experimentally 
observed one (7.13), if the polymer chain is assumed to bear an effective charge Z* = 
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(fN) 1 ' 2 . The particle concentration is replaced by c/N, and the hydrodynamic radius is 
proportional to N. The physical origin of the effective charge Z*, however, is unclear. 

This discussion indicates that a full theoretical understanding of the viscosity of poly- 
electrolyte solutions is far from being achieved. The main results can be summarized 
as follows. In all models, a maximum in rj r is predicted when the salt and the coun- 
terion concentration are of the same order of magnitude, fc ~ n. The maximum can 
occur, depending on chain length and salt concentrations, either in the dilute or in the 
semi-dilute region. For short chains or small salt concentrations, it thus takes place in 
dilute solution. The appropriate theory to describe the peak should be similar to that of 
references [109,108]. At concentration higher than the peak concentration, the reduced 
viscosity behaves as fNc^ 1 / 2 . Although the behaviour is similar above and below the 
overlap concentration, the interpretations that have been proposed are very different. For 
longer chains, or higher salt concentrations, fc ~ n corresponds to a semidilute solution. 
In that case, the viscosity can be described using equations (7.15) (for c larger than c*) or 
(7.21) (for c near c*). Below c*, the appropriate description is again that of [109,108]. 

8. Conclusions 

In this review we have presented the recent theories that describe the conformational 
and dynamical properties of charged polymer chains. Much of the focus has been put on 
the long range character of the electrostatic interactions and on the role of the small coun- 
terions that insure the electrical neutrality. Little attention has been paid to interactions 
of non electrostatic origin [110] which have only been introduced at some places in terms of 
virial coefficients. Most experiments on polyelectrolyte solutions are made in water that is 
not a good solvent for polyelectrolytes which are organic polymers. The non electrostatic 
interactions between molecules dissolved in water are not simple Van der Waals interac- 
tions and have components due to hydrogen bonding or to the hydrophobic effect [23] . It is 
not obvious that these very specific interactions can be modelled simply in terms of virial 
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coefficients. Experimentally these strong attractive interactions lead in some cases to the 
formation of aggregates which can be responsible for the slow modes [111,96,112] some- 
times observed when measuring the relaxation of concentration fluctuations via quasielastic 
light scattering. The simple theoretical models that we have described do not take into 
account these complicated effects which are in many cases dominant. This makes a direct 
comparison between the theory and the experimental results rather difficult and we have 
not attempted here any quantitative comparison. 

Despite this difficulty, some experimental features seem to be rather general in poly- 
electrolyte solutions. The structure factor of a polyelectrolyte solution always shows a 
peak at low ionic strength. The position of the peak has been discussed insection 6. both 
in dilute and semidilute solutions. Most experiments are also consistent with a rodlike 
behavior of the polyelectrolyte chains at small length scales. The rigidity of the chain is 
then characterized by an electrostatic persistence length. The variation of the persistence 
length with ionic strength or with the polymer concentration is however still a matter of 
controversy as seen in section 4. and more complete theories are certainly needed. The 
effect of heterogeneities of the charge distribution or of non electrostatic interactions on 
the persistence length do not seem to have been discussed. A good agreement between 
theory and experiment seems to exist only for rigid polymers where the theory of Odijk 
or Skolnick and Fixman [113,114] gives a quantitative description when counterion con- 
densation is properly taken into account. Fuoss law and the variation of the viscosity of 
polyelectrolyte solutions with the concentration have also been extensively confirmed ex- 
perimentally. The theoretical results presented in section 7. remain at the level of scaling 
laws and are not based on very systematic arguments. It should be also noted that the 
derivation of Fuoss law in semidilute solution based on the Rouse model implicitly assumes 
that there is only one length scale and therefore that the chain persistence length is of the 
order of the correlation length £ as suggested by some of the theoretical models. 
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An alternative way to study the properties of polyelectrolyte solutions is numerical 
simulations [12]. All the parameters are well controlled in the simulations and they should 
allow to sort out the respective roles of the electrostatic and non electrostatic interactions 
and thus to provide quantitative tests of the theoretical models. Many simulations have 
been performed on isolated polyelectrolyte chains in a salt free solution. There is in general 
quantitative agreement between the results and the electrostatic blob model. The simula- 
tions in the presence of salt often performed by using the Debye-Hiickel potential are less 
conclusive. For rigid polyelectrolytes, the agreement with the Odijk-Skolnick and Fixman 
theory is good; for flexible weakly charged polyelectrolytes, there is no good separation of 
length scales between the screening length and the persistence length and a quantitative 
test of the theories is difficult. Recently, Kremer and Stevens [115] have performed simu- 
lations on semidilute polyelectrolyte solutions, taking into account explicitely the discrete 
counterions. The polymers are strongly charged flexible chains in the vicinity of Manning 
condensation threshold. This is a very heavy numerical work but it should allow a very de- 
tailed description of the influence of electrostatic interactions in polyelectrolyte solutions. 
The extension of these simulations to weakly charged or more rigid polyelectrolytes for 
which the analytical theories have been constructed would provide strong tests of these 
theories. 

The review has also been limited to simple polyelectrolyte chains comprising only 
monomers with the same charge and neutral monomers. There are several other polymeric 
sytems where electrostatic interactions play an important role. Polyampholytes [116,117] 
are polymers that carry charges of both signs. Due to the attractions between opposite 
charges, they are often insoluble in salt free water and can only be dissolved at a finite ionic 
strength. They also show an antipolyelectrolyte effect, the reduced viscosity increasing with 
ionic strength. Polysoaps are comblike polymers, often polyelectrolytes, where hydrophobic 
side chains are grafted on the backbone [1,118]. They have properties combining those of 
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polyelectrolyte chains and of small ionic surfactants and form for example intramolecular 
micelles. Another system of interest is mixtures between neutral and charged polymers 
[119,120] which have been studied both experimentally and theoretically. The addition of a 
charged polymers can increase the solubility of a neutral polymer in water. The formation 
of mesophases is also expected in these systems. 

Throughout this review we also have implicitly assumed that the solvent (essentially 
water) is a continuous dielectric medium with a uniform dielectric constant. At high 
polymer concentration the dielectric constant crosses over from that of the polar solvent 
to that of an organic medium (essentially polymer) [121]. If the dielectric constant is low, 
charges of opposite signs are not dissociated and the polymer has an ionomer behavior 
dominated by attractive interactions between dipoles. Even in a more dilute solution, 
it is not clear that locally, the dielectric constant has the solvent value; this could turn 
out important in Manning condensation theory for example. There seems to be little 
theoretical work on these issues. Also, the small ions have been considered as pointlike, all 
specific interactions due to hydration or to the finite size of the ions have been neglected. 
Small differences between ions due to these effects are usually observed. 

One of the essential conclusions of this review is that polyelectrolytes remain a poorly 
understood state of matter and that much work is still needed both from the experimental 
and from the theoretical point of view to reach a degree of understanding equivalent to 
that of neutral polymers. ^From a theoretical viewpoint, polyelectrolytes are in most cases 
strongly coupled systems, involving many different length scales. It might well be that in 
such situations, an appropriate description of the system can only be obtained by the use 
of sophisticated liquid state theories, such as those recently developped for dense neutral 
polymers [122]. Experimentally, it is surprising that relatively few systematic studies on 
charged polymers, using modern investigation tools, in particular neutron scattering, are 
available. Such studies were initiated in the late seventies (see e.g [123]), but have not been 
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as conclusive as the equivalent studies on neutral polymers. Finally, it is obvious that the 
area of polyelectrolyte solutions ofters a number of problems that can be approached by 
simulation. Simulations that explicitly include counterions can be used to assess the accu- 
racy of the Debye-Hiickel or Poisson-Boltzmann approximations [124]. Density functional 
based simulation methods [125] for classical charged systems can be used to go beyond 
such approximations, even in the presence of added salt. Finally, the dynamical properties 
of charged polymers have not, up to now, been investigated in numerical studies. 
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APPENDIX A 
Effective interaction between charged monomers 

In this appendix, it is shown that the effective interaction between charged monomers 
can, within the linear response approximation, rigorously be written as a sum of pair- 
wise additive interactions (equation (2.16). As a byproduct of the calculation, the various 
stucture factors (ion- ion, polymer-polymer and polymer-ion) are linked within this approx- 
imation, by simple relations. 

The system is a polymer solution with a concentration c of (positively) charged monomers, 
and a concentration n of monovalent added salt. For simplicity we assume that the coute- 
rions are identical to the salt cations. The densities of charged monomers, charged positive 
and negative ions are denoted by pp(r), n+(r) and n_(r), respectively. By integrating out 
the coordinates of the small ions, the partition function of the system can formally be 
written as 

Z = Tr p exp(-(tf P + F t [i;])/k B T (A.l) 

where Tr p denotes an integration over monomer coordinates, and Hp is the energy of the 
polymer (2.1). Fi[ip] is the free energy of the small ions in a potential that is the sum of 
the external potential and of the electrostatic potential i/>(r) created by the polymer 

V(r) = k B Tt B J . (A.2) 

The free energy Fi is obtained by minimizing with respect to the ionic densities n+ (r) and 
n-(r) a free energy functional, which in the mean field or Poisson-Boltzmann approxima- 
tion reads 

k B T J dr{n + (r) ln(n+(r)) + ra_(r) ln(ra_(r))) 
+ \t B K B T J ^ (" + W-n4r)y)-n-(0) + J 

(A.3) 
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The linear response, or Debye-Hiickel, approximation, involves a further simplification of 
equation (A. 3), by expanding the first two terms to second order in the deviations of n+ 
and n- from their average values. The minimisation is easily carried out, and the resulting 
expression for Fj is 

with k 2 = 4n£ B (2c+n). The interaction that results from adding this contribution from the 
ionic free energy to the electrostatic repulsion between the monomers, J d k^ k i B pp(k)pp( — k), 
gives the screened pair potential (2.16). 

The Debye-Hiickel calculation can easily be extended to compute the response of vari- 
ous densities to external fields, or equivalently the partial structure factors. The results are 
simple in the case of a salt concentration much larger than the counterion concentration 
(n » c). Defining the ionic charge density as pz = n+ — n_, the following relationships 
are obtained between the Fourier transforms of the polymer-polymer , polymer-charge and 
charge-charge correlation functions (denoted by <Spp, Spz and Szz> respectively) 



Spz(k) = 
Szz(k) = 




(A.5) 



The total charge structure factor, Spp + Szz ~ ^Spz-, vanishes as 2nk ■ / k in the small k 
limit, as implied by the Stillinger and Lovett sum rules [22]. 
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APPENDIX B 
Relaxation and electrophoretic effects 

In this appendix, we derive equations (7.1) and (7.3), describing respectively the relax- 
ation and electrophoretic effects. As usual, the two effects are considered independently 
[91], since the coupling between them is relevant only at a level that goes beyond the limits 
of the Debye-Hiickel linearized theory. 

B.l Relaxation field 

The relaxation field is the electric field created by a charged particle (charge Qe) 
moving in a monovalent salt solution (concentration n) with a constant velocity V. As 
explained in section 7.1, this relaxation field arises from the deformation of the Debye- 
Hiickel polarisation cloud when the charge is set into motion. At equilibrium, the densities 
77.+ (r) of positive and n^r) of negative ions that surround the charged particle (placed at 
the origin) are 

n { ± q \r)=n(lTQyexp(- K r)) . (B.l) 

In a time dependent situation, however, the Boltzmann statistics that is used to obtain 
these densities must be replaced with a diffusion equation 

(0 dn± ^ t] = V • (fc B TVn±(r, t) =F n±(r, f )E(r, t)) . (B.2) 

where Co is the ionic mobility. The electric field E is the sum of the field Ep(r — Vt) 
created by the moving particle and of the electric field created by the ionic charge density, 
E^, which verifies Poisson equation 

V ■E i (r,t) = 47r£ B (n+(r,t)-n-(r,t)) . (B.3) 

The solution of (B.2) in the stationnary state is of the form n±(r,t) = n± q \r — Vt) + 
5n±(r — Vt). Inserting this form into (B.2) and (B.3), and linearizing with respect to V 
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and Q, the following equation for 5n+ — 8n- is obtained. 



(V 2 + ^ 2 )(^ + (r) - Sn-(r)) = -fi-V ■ V 

k B T 



£ 

(2nQ— exp(-Kr) 
r 



(B.4) 



The resulting ionic charge density can be expressed as 

( fa+ ( r )- fa _(r)(r) = ^/A (^^^expHk-r) (B.5) 
and the resulting electric field is given by equation (7.1). 
B.2 Electrophoretic effect 

We now consider the calculation of the hydrodynamic velocity field around a charge 
Q at the origin, when the charge experiences a force Q~E e xt created by an external electric 
field ~E e xt- The solvent experiences a force Q~E ex t5(r) created by the charge at the origin, 
and a body force e(n+(r) — n-(r))~E ex t. The equation that determines the hydrodynamic 
velocity field U(r) in the solvent is 

- Vs V 2 U + VP = QE ext 5(r) + e(n+(r) - n-(r))E ext , (B.6) 

where rj s is the viscosity and P the pressure. (n+(r) — n_(r)) is the charge density within 
the equilibrium Debye-Hiickel polarisation cloud. In Fourier space, the velocity field is 
given by a generalization of Oseen's formula, 

0i(k) = -^—^ (Sij - kikj/k^QBj). (B.7) 

Which gives in real space, 

U(r) = H s (r)(QB ext ) (B.8) 
with the screened Oseen tensor H s being 

^'H^T^ + ^il^- (B ' 9) 
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FIGURE CAPTIONS 



figure 1: schematic representation of the blob structure of a charged polymer chain. The 
chain is gaussian up to the scale stretched on larger length scales. 

figure 2: schematic "state diagram" for a polyelectrolyte solution, in the plane c (polymer 
concentration), n (concentration of monovalent added salt). The boundaries between the 
different regions of the plane indicate crossovers, not sharp transitions. Region I (dilute) 
is discussed in sections 2, 3 and 6.1. Regions II, III and IV (semidilute) are discussed in 
section 6.2 and 6.3. Region V (concentrated) is discussed in section 6.5. 

figure 3: Phase diagram of a dense polyelectrolyte solution in a bad solvent. Upper part: 
phase diagram in the concentration, temperature plane for a fixed salt concentration. 
Lower part: phase diagram as a function of salt concentration and temperature for a given 
polymer concentration. 

figure 4: schematic representation of the reduced viscosity as a function of polymer 
concentration c at a fixed salt concentration n. 
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